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ABSTRACT 


Observation of Bow Shock Protons at the Lunar Orbit 
by 

John Louis Benson 

Protons with energies ranging from about 500 eV to 
3500 eV are observed by the Suprathermal Ion Detector Ex¬ 
periment on both the dusk and dawn sides of the magneto¬ 
sphere. On each lunation these particles appear as a rather 
continuous phenomenon for 3 to 5 days after crossing from 
the dawn-side magnetosheath into the solar wind and for 
about 2 days prior to entering the dusk-side magnetosheath. 

Transverse proton flows of this type have been observed 
by Asbridge et al. r [1968], Scarf et al. / [1970] and are be¬ 
lieved to be deviated solar wind particles that are acceler¬ 
ated at the earth's bow shock and emitted outwardly along 
interplanetary field lines at greater than solar wind 
energies. 

Data from the SIDE and from the Explorer 35 lunar orbit¬ 
ing magnetometer have been analyzed and these data indicate 
that the transverse ion flows observed by the SIDE in the pre 
and post bow shock crossing regions of the lunar orbit are due 
to these deviated solar wind particles. 

A computer model based on EXB drift trajectories for 
particles leaving the shock has been developed and synthetic 
particle data produced by this model are in good agreement 
with the observed data. 
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Chapter 1 

This thesis is concerned with the study of the plasma 
environment of the moon; particularly when the moon is outside 
the earth's magnetosphere and in the interplanetary medium. 
Much has been learned about this environment as a result of 
placing particles and fields experiments directly on the lun¬ 
ar surface during the six Apollo moon landings and the follow¬ 
ing work consists of the analysis of data returned to earth 
by one such experiment, the Suprathermal Ion Detector Experi¬ 
ment, otherwise known as SIDE. 

The plasma environment of the moon should not be consid¬ 
ered independently of the earth's magnetosphere and the in¬ 
terplanetary medium and a brief review of the basic ideas and 
observations concerning the complex interactions of the solar 
corpuscular radiation, (solar wind), and the near extra-mag- 
netospheric region with the earth moon system is in order 
before going on to a detailed discussion of the results of 
this analysis. 

1.1 General Survey 

The sun continuously ejects a tenuous supersonic plasma, 
known as the solar wind, that consists mostly of protons and 
electrons in equal amounts. The plasma is tenuous in that 
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typical number densities are on the order of a few particles 
per cubic centimeter and supersonic in that the bulk flow 
velocity of the plasma is greater than the Alfven speed in 
the medium. The Alfven speed is, 

V A =B A 47r P l 1 / 2 ~ 50 km/sec 1-1 

in the interplanetary medium and is the quantity analogous to 
the sound speed in a collision dominated medium. The average 
bulk flow velocity of the solar wind is about 400 km/sec 
giving an Alfvenic Mach number, V bf /V A of about 8, indicating 
the highly supersonic flow. 

The solar wind drags the sun's magnetic field lines 
along with it due to the essentially infinite conductivity 
in the plasma in such a manner that the ambient interplanetary 
magnetic field lines intersect the earth's orbit at approxi¬ 
mately a 45° angle to the earth-sun line. The field lines 
in figure 1-1 are shown at this "garden'hose angle." 

The radial flow of plasma from the sun appears to come 
from a direction about 5° to the west of the solar direction 
due to the earth's orbital velocity of about 30 km/sec. 

Upon reaching the orbit of the earth, the solar wind 
plasma is unable to penetrate the earth's dipolar field; it 
therefore compresses the field creating a cavity in the flow 
known as the magnetosphere. The solar wind flows around this 



Figure 1-1. View of the magnetosphere and the 
lunar orbit in the ecliptic plane. Typical magnet¬ 
ic field lines in the interplanetary medium and in 
the magnetosphere are shown. 



Figure 1-1 
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cavity in a fashion resembling hydrodynamic flow. The 
boundary of the cavity is called the magnetopause and the 
distance from the center of the earth along the earth-sun 
line to this boundary is on the order of 10 R £ (earth radii). 

The supersonic nature of the solar wind gives rise to a 
detached shock wave a few R g upstream of the magnetopause. 

This standing shock/ known as the bow shock, was predicted by 
Zhigulev and Romishevskii in 1959 and by Oxford in 1962. Its 
existence was confirmed by the IMP 1 satellite, Ness [1964]. 
Due to the excessively long collision mean free path of a 
solar wind particle, — 1 A.U., the shock is not collision 
dominated or even collision influenced. Nevertheless, the 
high conductivity in the plasma causes it to behave somewhat 
like a fluid and shock calculations using the "fluid dynamic 
analogy", Spreiter et al. [1966], work fairly well in pre¬ 
dicting the location of the shock. See figure 1-2. A few 
R g upstream of the shock there is only relatively cool, 

(T ~ 10 °K) , tenuous, (n 10/cm ), unshocked plasma with 

a weak interplanetary field, (B ~ 57 ), imbedded in it. 

For a summary of the properties of the quiet solar wind, see 
table 1-1. Closer to the shock, where interplanetary field 
lines intersect the shock surface, complex interactions of 
the solar wind and the shock take place and energetic parti¬ 
cles have been observed flowing upstream of the shock in this 
region by several investigators? Asbridge et al. [1968], 
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Figure 1-2. Measurements of the positions of the mag¬ 
netopause and bow shock from IMP 1,[after Ness, (1965)]. 


-30 — 
K-' a ) 



Figure 1-2 
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Table 1-1 

Observed Properties of the Quiet Solar Wind 


bulk flow velocity 

320 km/sec 

number density 

8 cm 8 

proton temperature 

4X10 4 °K 

proton thermal anisotropy ratio 

2 

electron temperature 

1 to 1.5X10 5 °I< 

magnetic field strength 

57 

proton flux 

2.4X10 8 cm 2 sec" 1 

kinetic energy flux 

-2 -1 

.22 erg cm sec 

kinetic energy density 

7X10 9 erg cm 8 

proton thermal energy density 

6X10 erg cm 8 

electron thermal energy density 

1.5X10 10 erg cm - ' 

magnetic field energy density 

10 10 erg cm ^ 
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Prank and Shope £1968], Anderson [1968]and Freeman £1972]. 
These energetic particles, specifically protons in the energy 
range of .5 to 3.5 keV are the subject of this thesis. A 
more detailed discussion of these observations is given below. 
After passing through the shock, the temperature of the 

plasma is increased, {T ~10^°K), the magnetic field is com- 

-3 

pressed, (B ~ 20 y ), density goes up, (n ~ 30/cm ) and the 

bulk flow velocity goes down. Some of the ordered streaming 
energy is converted to disordered thermal energy, Figure 1-3 
summarizes the changes in plasma parameters across the shock. 
The region of space between the magnetopause and the bow 
shock is known as the transition region or the magnetosheath 
and the detailed plasma properties of this region near the 
lunar orbit are discussed by Fenner [1971,1974]. 

The inside of the magnetosphere is essentially shielded 
from the solar wind flow and figure 1-4 depicts the general 
particle and field morphology relevant to the Lunar environ¬ 
ment in this region of space. The earth's dipolar field is 
compressed on the day side in to about 10 R £ and drawn out on 
the night side at least several tens of R £ past the lunar 
orbit [Ness et al., 1967], in the anti-solar direction. 

The dipolar field lines are closed near the earth, but 
the distance in the anti-solar direction to which this clos¬ 
ing extends is not well established. At a distance of more 

than about 10 R , the field tends to point either toward or 
e 
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Figure l-3a. Typical changes in plasma parameters 
across the bow shock. Note that C=sound speed and 
b=Alfven speed. M Q =mach number =Alfven mach number, 
[after McKenzie, (1970)]. 

Figure l-3b. Measured changes in plasma parameters 
across the shock, [after Wolfe and McKibbin, (1968) 3 . 

Figure 1-4. View of magnetosphere in noon-midnight 
meridian plane, [after Ness, (1965)]. 
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away from the sun. Although the field lines may be weakly 
closed (B„ << 1) at 60 R e , the morphology of the geomag- 

netic tail in this region can generally be described as a 
planar neutral sheet with the field lines roughly parallel to 
the sheet above and below with the sense of the field being 
toward the sun north of the neutral sheet and away from the 
sun south of the neutral sheet. The thickness of the neutral 
sheet at 60 is not well known, but Bowling, [1973] has set 
a value of about 3 R £ at a distance of 30 R e in the anti¬ 
solar direction. 

The magnetic neutral sheet region is also filled with a 
relatively warm, {T^ ~ 10®°K), rarefied, {n ~ 10 /cm ), 
plasma that surrounds the neutral plane with a thickness of 
about 3 to 7 R e at the lunar distance as determined by CPLEE, 
[Rich, 1973]. See chapter 2 for a description of the Charged 
Particle Lunar Environment Experiment. The source of the par¬ 
ticles in this "plasma sheet" is not well established, but it 
is believed that they are solar wind particles that have 
entered the magnetosphere through the magnetopause or through 
open tail field lines at large distances in the anti-solar 
direction. 

The above brief description of the regions of the magnet¬ 
osphere and interplanetary medium through which the moon pass¬ 
es allows one to summarize the changing plasma environment of 
the moon, disregarding any lunar atmospheric phenomena, as 
follows. \ 
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Beginning at new moon, (see figure 1-1), the moon is in 
the undisturbed solar wind where the ordered streaming plasma 
rams the sunward side directly without going through a shock 
wave and the earthward side faces a plasma void that lines up 
with the solar wind aberration direction. The absence of a 
lunar bow shock was discovered from Explorer 35 magnetometer 
data, (Colburn, 1967] and the plasma void was identified using 
a plasma probe flown on the same spacecraft, [Lyon et al., 

1967] . The Explorer 35 magnetometer has also determined that 
a slight enhancement of the interplanetary field occurs while 
the satellite is in the moon's optical shadow, [Colburn et al., 

1968] . 

At first quarter the conditions are essentially the same 
as at new moon, as long as the interplanetary field does not 
change its streaming angle appreciably. That is, at first 
quarter the interplanetary magnetic field lines still do not 
connect the moon and the bow shock so that the bow shock can 
not yet influence the lunar environment. 

When the moon proceeds toward the dusk magnetosheath, it 
will begin to encounter field lines connected to the bow shock 
with the probability of observing bow shock associated part¬ 
icles becoming larger as the moon moves closer to the shock. 

After passing into the dusk magnetosheath, the plasma 
environment changes considerably from that just outside the 
shock. As discussed above, the density increases, the bulk 
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flow velocity decreases and the temperature increases. The 
magnetic field is enhanced in this region and particle spec¬ 
tra become very stable with peak energies around 500 eV. 

By the time the moon has reached full phase it is well 
into the geomagnetic tail where it experiences multiple 
.neutral sheet crossings and encounters with the plasma sheet; 
densities decrease from magnetosheath values and temperatures 
increase. 

As the moon passes from the full to the last quarter 
phase it crosses the dawn magnetosheath wheia the plasma par¬ 
ameters are similar to those in the dusk magnetosheath but 
differences in turbulent effects are observed. The dusk 
magnetosheath exhibits a correlation in turbulence and geo¬ 
magnetic activity (Kp) and the dawn magnetosheath does not. 
[Fenner, 1974]. 

The region between the dawn bow shock and the point at 
which the last interplanetary field line tangent to the shock 
intersects the lunar orbit, (see figure 1-1), is the primary 
region of interest in this thesis since the probability of 
observing bow shock protons is greatest here. 

1.2 Observations of bow shock associated phenomena, up stream 
of the shock . 

Several investigators have observed energetic particles 
and wave disturbances upstream of the earth's bow shock. A 
brief summary of those observations follows. 
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Waves of frequency 10-100 Hz propagating along 
interplanetary field lines have been detected upstrea of the 
earth's bow shock by the 0G0-1 search coil magnetometer, 

[Smith et al., 1967]. These waves are believed to be precur¬ 
sor waves propagating upstream in the whistler mode parallel 
to the field lines. Similar observations of whistler waves 
upstream of the shock were made by Heppner et al. [1967] 
although the frequencies were about 1 Hz with the waves being 
detected in small packets of 4 to 6 cycles. The amplitudes 
of these waves were found to diminish with increasing distance 
from the shock which led Heppner to suggest that they were 
generated at or were associated with the shock itself. 

A comprehensive study of upstream waves was made by 
Fairfield [1969], using Explorer 34 data. He determined that 
waves of frequency 0.01 to 0.05 Hz upstream of the bow shock 
are a permanent feature of the region of space extending 
from the bow shock out to about 34 R g . These waves are 
observed 90% of the time when the interplanetary field inter¬ 
sects the shock and Fairfield believes them to be locally 
generated Alfven waves due to the outward streaming protons 
observed by Asbridge and other investigators. He argues that 
the whistler waves must have a fairly high frequency 
(v ^ 10 Hz) to propagate upstream as shown by Smith et al. 
[1967]. Processes that would explain the necessarily large 
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Doppler shifting to .01 Hz in the rest frame of the spacecraft 
would be difficult to explain. 

Greenstadt et al. [1968] have contributed further to the 
association of these low frequency waves with particles by 
showing a correlation in magnetic field oscillation and 
particle fluxes outside the shock. Energetic electrons 
(energy ^40 keV) have also been observed in the upstream re¬ 
gion by Anderson [1968] , and by Jokipii [1968] . Protons of 
1 keV were observed upstream of the shock by Frank and Shope 
[1968 ] r with the IMP-4 satellite while the detector was look¬ 
ing in the anti-solar direction. Freeman [1972] has also 
observed energetic ions in the deep lunar night with the 
suprathermal ion detector on the lunar surface when the detec¬ 
tor was looking as near as 17° to the anti-solar direction. 
Freeman points out that these ions are probably not from the 
earth's bow shock and that they are most likely associated 
with deviated solar wind flow around the moon. Examples of 
outwardly flowing protons are reported by Asbridge et al. 
[1968 ] r using data from six Vela satellites. The results 
reported in their paper are particularly relevent to the 
present work and this paper will now be discussed in more de¬ 
tail. 

Asbridge and co-workers report the results of the 
analysis of particle data from the six Vela satellites which 
are in nearly circular earth orbits at 18 R g . The energy 
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analyzers on board are capable of detecting protons in the 
range of 0.1 to 20 keV. All of the Vela satellites are spin 
stabilized and sample only certain directions in spacecraft 
azimuth. Vela's 2 and 3 spin at about 2 rps and take five 
measurements near the solar direction and others at 90°, 180% 
and 270° from the solar direction. The Vela 4 satellites spin 
at about 1 rpm taking seven measurements of proton spectra 
near the solar direction and a measurement each 22 1/2° of 
rotation thereafter. 

Several examples of proton flows in the energy range 3 to 
6 keV are presented in this paper. These protons are observed 
to come from the direction of the earth's bow shock and figure 
l-5a,b,c shows two examples of the dawn side flows and one 
example of a dusk side flow. An analysis of four months of 
Vela 3A and 3B data shows that these outwardly flowing protons 
are a permanent feature of the dusk-dawn bow shock region 
and figure l-5d summarizes these results. Arrows indicate 
direction of flow and dots indicate flow in the solar direc¬ 
tion. A higher incidence of outwardly flowing protons is 
noted when the spacecraft are on the dawn side .of the magnet¬ 
osphere. The authors suggest that these protons are solar 
wind particles that are accelerated at the bow shock and 
emitted outwardly from the shock along interplanetary field 
lines. Figure 1-6 summarizes the model suggested by these 
investigators and it is shown in this thesis that this model 
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Figure l-5a. Dawn-side transverse proton flow from 
the bow shock as observed by Vela 3B; 291° spacecraft 
azimuth points toward the bow shock while on the dawn 
side, (energy 4 keV). 

Figure l-5b. Dusk-side transverse flow from Vela 2B 
observations at 18 R e> (energy ax 3 keV). 

Figure l-5c. Vela 3A dusk-side observation of transverse 
flow at 89° front the solar direction, (energy ax 5 keV) . 

Figure l-5d. Summary of four months of Vela 3A and 3B 
data. Arrows indicate direction of the flow; dots indi¬ 
cate flow in the solar direction. 

Figure 1-6. Diagram of Asbridge model to explain the 
observed transverse proton flows. 

.Figures 1-5 and 1-6 are from [Asbridge et al. (1968)]. 
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Figure 1-6 
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is at least partially correct. A more realistic model is de¬ 
veloped and discussed along with the Asbridge model in chap¬ 
ter 3. 


1.3 Research Objectives 

The SIDE observes a variety of ion spectra in the region 
of space between the bow shock and the terminator. The dom¬ 
inant spectrum is one that is depleted in the lower energies 
and peaks toward the higher energies. Other spectra show 
double peaks with one doubly peaked type occuring near the 
bow shock and the other near the terminator. The lack of a 
complete understanding of the ion data in this region and the 
fact that the dominant spectrum is what one would expect from 
bow shock protons lead to the following research objectives; 

a. ) to determine whether or not the moderate to higher 
energy ions observed by SIDE for several days after passing 
out of the magnetosheath and into the solar wind are parti¬ 
cles from the earth's bow shock. 

b. ) to identify these particles in the data, further 
explaining some of the more complicated spectra observed in 
this region , [Lindeman, 1971]. 

c. ) to propose and develop a theoretical model to ex¬ 
plain how bow shock protons can reach the lunar orbit. 
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Chapter 2 

2.1 Basic Theory of instrument Operation 

.e basic idea behind the SIDE and most particle detectors 
with energy discrimination capability is simply that the traj¬ 
ectory of a charged particle is altered upon its entering a 
region in which there exists an electric field. Detectors 
that are based on this principle often employ a set of parall¬ 
el plates across which a potential difference is applied to 
generate the electrostatic field. A channel electron multi¬ 
plier is placed at the end of the plates as a counter. If the 
configuration of the plates is such that the electric field is 
applied perpendicular to the trajectory of the incoming part¬ 
icle, then the amount of deviation from that initial trajec¬ 
tory depends on the energy of the particle and the strength of 
the electric field. 

The very simplest configuration for such a detector would 
be a pair of straight parallel plates of length L that are 
separated by a distance d. Then if a potential V is applied 
across the plates, a particle entering the region between the 
plates on a trajectory that is parallel to the plates and at 
a distance of d/2 above the bottom plate, must have an energy 
that is at least equal to 


£= q L?V 2-1 

2d 2 - 

to get through the plates, where q is equal to the charge of 
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the particle. 

This configuration is not very useful since it only pro¬ 
vides the lower limit on the particle energy, i.e. a particle 
with energy greater than can also make it through the plates. 
The energies of particles that get into the counter can be 
confined to a narrow range if the parallel plates are curved. 

The SIDE employs such a pair of curved plates in both the 
total ion detector (TID) and the mass analyzer detector (M&). 

The potential, as a function of r, between the curved 
plates is given by, 

dbfr) = V ln(r/a) 2-2 

In (b/a) 

where b is the radius of curvature of the outer plate and a is 
the radius of curvature of the inner plate{see figure 2-1), 

V is the total potential difference across the plates. 

The electric field is given by, 

E = -V<p 2-3 

so that the field strength is 

E = V i 2-4 

In(b/a) r 

The force exerted on a particle of charge q, moving 
through the curved plates is , 
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Figure 2-1. Schematic representation of parallel 
plate electrostatic energy analyzers, flat and curved. 
Crossed E and B fields (Wien velocity filter) preceed 
the Low Energy Curved Plate Analyzer to form the Mass 
Analyzer. 
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and the centrifugal force on the particle of mass m, is 


If these forces are balanced, the particle will make it 
through the plates and into the counter. That is, the energy 
per unit charge of a particle that gets through is a function 
of the potential difference across the plates and the plate 
separation. 

C, = - 2-7 

2 In (b/a) 

Equation 2-7 is only approximately correct since equation 
2-2 was derived for infinitely long cylindrically concentric 
plates. The parallel plates of the detector are merely a sec¬ 
tion of the concentric cylinders and hence, edge effects are 
important in determining the correct equation relating the en¬ 
ergy and the potential difference across the plates. The geo¬ 
metry of the plates and these edge effects give rise to a fin¬ 
ite range of energies that may be passed at a particular val¬ 
ue of V. The results of the calibration of the Apollo 14 in¬ 
strument to determine the size of the energy passband are 
shown in figure 2-2. The bandwidth at the nominal center 
energy is approximately 10% of the center energy. 

2.2 Mass Analyzer 

The curved plates of the mass analyzer are preceeded by a 
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Figure 2-2. Some calibration results of the Apollo 14 
Total Ion Detector. Upper part of the graph shows the 
unidirectional geometric factor as a function of the 
nominal center energy. Lower part shows the energy band¬ 
pass (FWHM) as a function of the nominal center energy. 
Note that the bandpass is approximately 10% of the center 
energy. 
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Wien velocity filter which consists of crossed electric and 
magnetic fields that select the incoming particle’s velocity 
before it is energy-analyzed in the curved plates. If the par¬ 
ticle’s velocity and its energy per unit charge are known, then 
the mass per unit charge is simply, 

m/q = 2 £*/(v 2 q) 2-8 

Upon entering crossed electric and magnetic fields, a 
charged particle experiences two forces, one associated with 
the E field and one associated with the B field. The relative 
strengths of these two fields determines the resultant trajec¬ 
tory of the particle. The force due to E is, 

P e = q E 2-9 

and that due to B is, 


so that if, 


2-10 

~ F e =0 2-11 

the particle goes through the fields undeviated from its ini¬ 
tial trajectory. This requires that the particle have a velo¬ 
city, 

v = E/B 2-12 


to pass through the velocity filter and into the curved plate 
analyzer. 
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The E field of the velocity filter is provided by flat 
parallel plates and the B field is provided by a permanent mag¬ 
net. By adjusting the voltage on the plates, various ratios 
of*E to"B can be obtained, thus allowing only particles of cer¬ 
tain velocities through the filter. 

The basic idea behind the mass analyzer is also depicted 
in figure 2-1. 

2.3 Description of .the instrument 

Figure 2-3 shows schematically the arrangement of the var¬ 
ious components of the instrument. The SIDE has two positive 
ion detectors, one of which is preceeded by a Wien velocity 
filter, as discussed above, and this detector is known as the 
mass analyzer detector (MA). The curved plate analyzer, asso¬ 
ciated with the mass analyzer is connected to a voltage step¬ 
ping supply and is sensitive to energies per unit charge rang¬ 
ing from 48.6 to 0.2 eV/q and is referred to as the Low Energy 
Curved Plate Analyzer (LECPA). After passing through the 
velocity filter and the LECPA, an ion enters a counter known 
as a Bendix Channeltro®. This channel electron multiplier is 
held at a potential of -3.5 kV to accelerate the ions after 
they pass out of the LECPA. This post acceleration of the ions 
affords a higher detection efficiency of the channeltron which 
emits a cascade of secondary electrons whenever an ion strikes 


its interior. 




Figure 2-3. Schematic diagram of the electrical con¬ 
nections of the TID and MA. [Hills et al. (1971)]. 






Figure 2-3 


LUNAR SURFACE 
GROUND PLANE 









25 


The curved plate analyzer associated with the total ion 
detector (TID), is known as the High Energy Curved plate Analy¬ 
zer (HECPA). It is also connected to a voltage stepping sup¬ 
ply and is designed to pass energies per unit charge ranging 
from 3.5 keV/q to 10 eV/q. Ions leaving the HECPA are also 
post accelerated into a channeltron. The physical arrangement 
of the various components of the instrument is shown in figure 
2-4. 

The velocity filter for each of the three instruments 
is different and hence each measures a different mass range. 
Table 2-1 summarizes the energy and mass sensitivities of the 
three instruments and figure 2-5 shows some results of the 
mass analyzer calibration. 

Each SIDE is also equipped with a spider web-like screen 
(see figures 2-6,2-7) made of conducting material, that con¬ 
nects the instrument bo the lunar surface through a voltage 
stepping supply. Figure 2-7 is a schematic representation <f 
this connection. The voltage supply is also connected to a 
wire grid above the ion entrance apertures. This grid is 
at instrument ground. 

The purpose of the ground plane stepper is to compensate 
for any existing lunar surface potential and to determine its 
value. For example, if the lunar surface potential is +10 
Volts with 0 volts across the stepping supply and the instru¬ 
ment is electrically connected to the surface, then a 10 eV 
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Figure 2-4. Diagram of the internal configuration of 
the various components of the SIDE, [after Freeman et al. 
(1969)]. 

Figure 2-5. Results of mass analyzer calibration, 
(plotted from data compiled by Lindeman [1973]). 

Figure 2-6. External view of SIDE showing ground plane ' 
and grid over the entrance aperature of the instrument, 
[from Freeman et al. (1969)]. 

Figure 2-7. Schematic diagram showing relation between 
grid and ground plane. The ground plane is tied to the 
aperature grid through a stepping voltage supply. 
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proton would not be detected due to the +10 V potential at 
the entrance aperture. If, on the other hand, the voltage 
supply maintained a -10 V potential on the grid, with respect 
to the ground plane, the net potential seen by an incoming 
10 eV proton would be zero and it would enter the detector as a 
10 eV particle and would be counted as such. This connection 
between the lunar surface and the instrument through the step¬ 
ping voltage supply can also be used to determine the value 
of the lunar surface potential. If we assume that the ambient 
lunar atmosphere consists of essentially thermal ions (energy 
0 eV) then as the ground plane is stepped through various 
voltages, the ions will begin to appear in the energy channel 
corresponding to the absolute difference between the lunar 
surface potential and the voltage on the grid. For example, 
if the voltage on the grid is -2C V and particles show up in 
the 10 eV/q channel, then we can expect that if these particles 
were indeed from the ambient lunar atmosphere, having almost 
no energy initially, that the lunar surface potential is about 
+10 V. 

The usefulness of the ground plane stepper in determining 
the lunar surface potential has been demonstrated by Fenner 
et al. [1973], and it has been determined that a potential of 
about +10 V is appropriate for the dayside of the moon. Based 
on data collected near the terminators a value of about -100 V 
has been set for the night side. Figure 2-8 shows the 
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Figure 2-8. Results of determination of the lunar 
surface potential in two regions of the orbit. 
[Fenner et al. (1973)]. 
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lunar surface potential data thus far obtained. Theoretical 
predictions of this potential, based on various values of 
photoelectric emissivity of the lunar surface, have been made 
by Manka [1972] and these theoretical predictions are also 
displayed in figure 2-8. 


2.4 Flux measurement 

The counting rate R, which is the quantity actually avail¬ 
able from SIDE data, can be converted into a more useful quan¬ 
tity, namely the unidirectional differential energy flux 

j 

The counting rate at energy £ is given by, 


R = ( fg< € 

where, 


e , 0 ) j ( £ , Q , <P )&6 d<£ d£ 


2-13 


g ( € r 6 , <p) =unidirectional geometric 
factor at energy £ 


and 


j ( £ , 8 , <P ) =unidirectional differential 


flux 

If we assume that the flux is independent of 6 and <p over 
the solid angle of the detector's field of view so that 
j { £ r 9 r <f> ) is approximately constant, then. 


~ I 3 ( € * 6 , <p) I g { £ , 9 , <p) dd dtp df 
= /j(£,Q)G(£)df 


2-14 
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where, 

G(f)=/g(f,0,0)<aa=E(£)e(f) 2-15 

and H C £)=*(£) 0(£) and e ( £ )=counter efficiency. 
eff eff 

a is the effective area of the entrance aperture as seen 
eff 

by an incoming particle and fi e j=£ is th e effective solid angle. 
If we also assume that the flux is constant over the energy 
bandwidth of the detector, then 
R= j(fffl) /G(f)d£ 

J € 

R = j ( f , a ) G q 2-16 

G q = / G ( £ ) d£ (cm 2 ster ev) 

where the integration is carried out over the energy bandwidth. 

Therefore, the unidirectional differential flux at energy 
£ is given by. 


-2 -1 -1 “1 

j ( £ , n )=R/G q (cm sec ster eV ) 2-17 

likewise, the omnidirectional differential energy flux is, 
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Figure 2-2 shows the results of the geometric factor cal¬ 
ibration for the Apollor XXV instrument [Lindeman, 1973] . 


2.5 Deployment 

The SIDE is part of a group of experiments, known as the 
Apollor Lunar Surface Experiments Package (ALSEP), designed 
to investigate various phenomena associated with the moon it¬ 
self and its interaction with the solar wind and the earth's 
magnetosphere. SIDE'S have flown on three Apollo missions; 
all were successfully deployed and all are still operating 
at the time of this writing. 

The Apollo XII site is located at about 3°S latitude 
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and 23°W longitude in Oceanus Procelarum (Sea of Storms). 

The Apollo 12 SIDE was placed on the lunar surface by Alan 
Bean and Charles Conrad during the first EVA which began at 
about 6 a.m. EST November 19, 1969. 

The Apollo XIV mission took place in early 1971, result¬ 
ing in the placement of the Apollo 14 SIDE on the lunar sur¬ 
face at about 3°S latitude and 17°W longitude in the Fra Mauro 
formation. Alan Shepard and Ed Mitchell deployed the ALSEP 
on their first EVA which began about 9 a.m. EST, February 5, 
1971. 

The Apollo XV SIDE was deployed by Dave Scott and James 
Irwin in the Hadley-Apennine region of the moon at about 26 °N 
latitude and 3°E longitude again during the first EVA which 
began about 8 a.m. EST July 31, 1971. 

Figure 2-9 shows the positions of the lunar surface of 
the landing sites and their seleonographic coordinates. 
Photographs of the three instruments are shown in figure 2-10. 

The Apollo 12 and 14 SIDE'S were equipped with a bubble 
leveling device attached to the top of the instrument that 
allowed the astronauts to align the vertical axis of the SIDE 
with the local vertical to within a 5° accuracy so that the 
look direction of the detector is approximately in the eclip¬ 
tic plane. The Apollo 15 instrument was equipped with an ad¬ 
justable leg that allowed compensation for the high latitude 
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Figure 2-9. Seleonographic locations of Apollo 12, 

14 and 15 ALSEP sites. 

Figure 2-10. From left clockwise, photographs of 
Apollo 12 and 14 SIDE'S deployed on the lunar surface 
Apollo 15 SIDE in the lab. [after Medrano, 1973]. 





MOON 


NORTH 



SOUTH 


LONG LAT M D Y 
APOLLO 12 23.4W 3.2S II 19 69 

APOLLO 14 17.4 W 3.6S 2 5 71 
APOLLO 15 3.6E 26.IN 7 31 71 


Figure 2-9 
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deployment so that its detector look direction is also approx¬ 
imately in the plane of the ecliptic, (see figure 2-lla) 

The look direction of the actual detector in each of the 
instruments is at a 15° angle from the instrument x-axis, which 
coincides with the local vertical at the ALSEP site for the 12 
and 14 instruments. This 15° offset determines the actual look 
direction, in longitude, with respect to the moon. The config¬ 
uration in longitude is illustrated in figure 2-llb. Figure 
2-12 shows the body axes of the three instruments. The relat¬ 
ionship of the body axes to the local coordinate system at 
ALSEP is as follows. 

The local coordinate system is known as system 1, and in 
the discussion in chapter 4, it is referred to as S'". It is 
defined by an x-axis that is aligned with the local vertical 
at ALSEP, a z-axis that is perpendicular to the x-axis and dir¬ 
ected northward along the meridian of deployment. The y-axis 
completes the orthogonal right handed system. That is, y 
points east. 

In the case of the Apollo 12 and 14 SIDE'S the body x-ax¬ 
is of the instrument coincides with the x-axis of system 1. 

The body x-axis of the 15 instrument is in the plane of the 
meridian of deployment but tilted southward by 26°. The body 
y-axis of the 15 instrument coincides with the y-axis of sys¬ 
tem 1 and the z-axis is orthogonal to x and y. The body z-ax- 
is of the 14 instrument coincides with system 1, but the body 
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Figure 2-11 


Figure 2-11, 


Figure 2-12. 


a. shows the configuration in latitude 
of Apollo's 14 and 15. Apollo 12 is ap¬ 
proximately represented by the drawing 
for Apollo 14. 

b. shows the configuration in longitude 

of all three SIDE’S. Apollo 14 SIDE looks 
most nearly toward earth. 

Relationship of the instrument body axes 
to the local ALSEP coordinate system 1. 




NORTH 
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A 

y- and z-axes of the 12 SIDE are the negatives of system 1. 

For further discussion of the relation of the body axes 
to other relevant coordinate systems, see chapter 4 and Hills 
[1969] . 

The deployment of the instruments in this fashion allows 
a directional sampling of magnetospheric and solar wind part¬ 
icles at any given time. The look directions of the three de¬ 
tectors for various positions in the lunar orbit are shown in 
figure 2-13. Note that the Apollo 14 SIDE looks almost direct¬ 
ly back at the earth at any position in the orbit, which makes 
this instrument particularly suitable for the observation of 
particles streaming outward from the earth's bow shock. Be¬ 
cause of the excessively high background counting rate in the 
Apollo 15 instrument and the discontinuous nature of the Apollo 
12 data for the period of interest in this thesis, only data 
from the Apollo 14 SIDE were used. 

2.6 Magnetometer [Colburn, 1969] 

The data on the interplanetary magnetic field and those 
from the geomagnetic field were collected by the NASA Ames mag¬ 
netometer onboard the Explorer 35 satellite which achieved lun¬ 
ar orbit on July 17, 1967 with a periselene of 1.4 R L (lunar 
radii), an aposelene of 5.4 and an orbital period of 11.5 
hours. 

The magnetometer contains three orthogonal fluxgate sen- 





Figure 2-13. Longitudinal look directions of the 
three SIDE’S for various positions in lunar orbit. 
Note that the Apollo 14 SIDE looks most nearly toward 
earth. Different sunset and sunrise times are due to 
the difference in longitude of deployment. 


[Hills et al. (1971)1. 
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sors which measure the ambient magnetic field in the neighbor¬ 
hood of the spacecraft and a vector measurement of the field 
is obtained every 6.14 seconds in the 20 gamma range. The 
signals from the sensors are spin-demodulated onboard to pro¬ 
vide a vector in a fixed spacecraft coordinate system, with 
the final data available from Ames in Solar Equatorial or Solar 
Magnetospheric coordinates, averaged over 81.8 second intervals 
The data are further processed on the SDS 910 computer at 
Rice to provide the field components in Solar Ecliptic coordi¬ 
nates , [Lindeman,1972]. 

The solar ecliptic coordinate system is defined as follows 
The x-axis is defined by the earth-sun line with positive x 
toward the sun. The z-axis is normal to the plane of the ec¬ 
liptic and positive z is toward north. North is the direction 
defined by the direction of the orbital angular momentem vec¬ 
tor of the earth, y completes the right-handed system. Solar 
ecliptic longitude is defined as 

0= tan" 1 (y/ x ) 2 " 21 

and solar ecliptic latitude is defined as 

6= tan -1 (z/(x 2 + y 2 ) 1/2 ) 2 ~ 22 

Due to a failure in the sun sensor in the spinning space¬ 
craft, the values of the solar ecliptic longitude are often 
in error. It is possible, however, to get a rough correction 
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for this error by observing the longitude data in the geomag¬ 
netic tail region where stable neutral sheet crossings are ob¬ 
served. 

If a neutral sheet crossing occurs when the spacecraft is 
near 180° in longitude, one would expect that while in the 
southern lobe of the field, the solar ecliptic longitude of 
the field would be about 180° and after the crossing, the 
field would have a longitude of approxiamtely 0°. By observ¬ 
ing the actual measured values and knowing the nominal values 
of the field longitude, the correction angle can be determined. 
An example of this is shown in figure 2-14 where a correction 
angle of about 28° was obtained for the lunation with which 
this thesis is concerned. For the present analyses, the un¬ 
corrected values of the longitude were used. Any physical 
implications of using these uncorrected longitudes are dis¬ 
cussed in chapter 5. It should also be noted that the Explor¬ 
er 35 magnetometer has no onboard recorder and therefore there 
are frequent data gaps of about 5 hours during each orbit due 
to the passage of the satellite behind the moon. 

2.7 Other lunar surface experiments 

The SIDE is only one of a variety of instruments design¬ 
ed to investigate plasma properties in the vicinity of the 
moon. These other lunar surface experiments provide data 
which can potentially be correlated with SIDE observations. 
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Figure 2-14. Solar ecliptic longitude of the magnetic 
field near the anti-solar position in lunar orbit. The 
technique for determining the correction angle is illus¬ 


trated. 





won 3S 
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Although detailed discussions of chese other experiments are 
beyond the scope of this thesis, the descriptions below are 
included to provide an i.nsight into the flavor of the particles 
and fields research being conducted on the lunar surface. 

2.8 CPL'E'F 1 -Charged Particle Lunar Environment Experiment 
The CPLEE consists of two charged particle analyzers, 

each containing five C-shaped channeltrons as negatively 
charged particle detectors and one channeltron as a positively 
charged particle detector, with electrostatic deflection plates 
providing the energy discrimination capability. Particles 
with energies from 50 eV to 50 keV can be detected by CPLEE 
from two directions simultaneously; one analyzer looks verti¬ 
cally with a 4°X20° field of view and the other detector is 
directed at 60° from the local vertical with the same field 
of view. 

The CPLEE has observed such phenomena as a photoelectron 
layer at the lunar surface, solar wind particles and penetra¬ 
tion of particles of auroral energies to lunar distances in 
the geomagnetic tail [O'Brien, 1971] . 

2.9 CCGE-Cold Cathode Gauge Experiment [Johnson et al., 1971] 
The cold cathode gauge experiment is designed to measure 

the neutral gas pressure of the lunar atmosphere. The instru¬ 
ment is attached to the SIDE and is connected to the ALSEP 
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central station through the SIDE, but it is a separate exper¬ 
iment consisting of a chamber in which a spool shaped cathode 
is connected to an electrometer. The cathode is surrounded 
by a cylindrical anode that is held at +4500 volts. When 
neutral gas enters the chamber, a self sustaining discharge 
occurs with the ion current through the cathode being a measure 
of the gas pressure. 

Results from the CCGE set neutral gas number densities 
5 —3 7 —3 

for lunar nightime of about 10 cm and about 10 cm for 
daytime with the daytime value partly due to local experiment 
outgassing. 

2.10 LSM--LPM-Lunar Surface Magnetometer and Lunar Portable 
Magnetometer. iDyal et al., 1972] 

The lunar surface magnetometers each consist of three 
orthogonal fluxgate sensors, mounted on booms, with the elec¬ 
tronics package located at the junction of the booms. This 
configuration provides a vector measurement of the static 
as well as the time varying magnetic fields at the ALSEP sites. 

The lunar portable magnetometers are self contained and 
also provide vector measurements of the field at several dif¬ 
ferent locations along the astronaut traverses. Table 2-2 
shows the range and sensitivity of the Apollo 16 LSM and LPM 
[Dyal et al., 1972]. Several of these magnetometers have 
flown on Apollo missions and a summary of the preliminary 
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Summary of Lunar Surface Magnetic Field Measurements 



Coordinates, deg 

Field magnitude, 7 

Magnetic field components, 7 

Sire 

n up 

B cast 

D north 

• Apollo If ' ? 

A LSI: P site 

8.9 S 15.5 E 

231 

*174 

-73 

134 

Site 2 


180 

-179 

-10 

19 

Site 5 


125 

120 

-25 

-25 

Situ 13 


313 

-141 

-209 

-186 

f inal Rover parking site 


121 

-42 

-93 

66 

Apollo 15 ALSEP site 

26.1 N 3.7 I: 

6 ±4 

4 +4 

1 +3 

4 ±3 

Apollo 14 landing site 

Site A 

3.7 S 17.S W 

103 ±5 

-93 ±4 

38 ±5 

-24 ±5 

Site (' 


43 ±6 

-15 ±4 

-36 ±5 

-19 ±8 

Apollo 12 Al-SI-P site 

3.2 S 23.4 W 

38 ±2 

-25.R i 1.0 

11.9 ±0.9 

-25.8 ± 0.4 


a-flie error in all components of the Apollo 16 field readmes is ±25 7. These limits will be reduced after luither data 
reduction. 


Table 2-3 
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Estimated 4 lie Flux Averages During the Times 
of Foil Exposure 


Mission 

| Exposure initiation | 

Exposure duration, 
hr: min 

A verage solar wind 

4 He flux. 

1 0 6 em'~scu'l 

Date 

(7.IU./., 
hr: min 

Apollo 11 

July 21, 1969 

03:35 

01:17 

6.2 1 1.2 

Apollo 12 

Nov. 19, 1969 

12:35 

18:42 ! 

8.1 J 1.0 

Apollo 14 

I'cb. 5, 1971 

15:15 

21:0(1 

4.2 i 0.8 

Apollo 15 

July 31, 1971 

19:36 

41:08 

17.7 i 2.5 • 

Apollo 16 

Apr. 21,1972 

23:01 

45:05 

12.0 t 1.8 


Solar Wind Abundance Ratios Determined From 
the Foil Exposure Periods of the Apollo Missions 


Mission 

4 He/ 3 He 

4 Hc/ 20 Nc 

20 Np /22 Ne 

22 Ne/ 21 Ne 

Apollo 11 

I860 ±140 

430 ±90 

13.S + 1.0 

_ 

Apollo 12 

2450±100 

620 ±70 

13.1 ±0.6 

26 ±12 

Apollo 14 

2230±140 

550 ±70 

13.65 ±0.50 

- 

Apollo 15 

2310±120 

550+50 

13.65 ±0.30 

31 ±4 

Apollo 16 

2260 ±100 

570 ±50 

13.80 ±0.40 

31 ±4 


Preliminary Results on Argon 
Abundances in the Solar Wind 


Mission 

20 Ne, 3 »Ar 

Apollo 14 a 

"VS 


Apollo 15 



Apollo 16 

29 ± 6 



First Results From the Analyses of the A luminum Foil 
From the Apollo 16 SWCExperiment 


Sample 

S’ 

4 He concentration, 
I0 J0 atoms/cm 2 

4 He/ 3 lle 

4 ile/ 20 Nc 

20 tXe/ 22 Ne 

22 Nc/ 2I Ne 

20 Ne/ 36 Ar 

Apollo 16: 








a 2-l 

5.02 

108 

2290 

510 

14.0 

32 

29 

2-2 

4.98 

107 

2300 

490 

13.9 

27 

30 

“2-3 

5.04 

101 

2280 

500 

13.5 

31 

27 

2-4 

5.17 

108 

2280 

490 

' 13.7 

27 

31 

“2-5 

10.35 

100 

2310 

490 

13.9 

35 

28 

Apollo 15: 








a 3-9 

3.13 

140 

2370 

480 

13.8 

33 

19 




Table 2-4 
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results of ambient field measurements at different locations 
is shown in table 2-3. 

2.11 S WCE-Solar Wind Composition Experiment iGeiss et al. f 
1972] 

The SWCE is one of the simpler lunar surface experiments 
and consists of exposing an aluminum foil strip, 30X140 cm, 
directly to the solar wind and then returning the foil to 
the laboratory for analysis. The foil is baked out in a 
vacuum chamber and in some cases melted to release the trapped 
solar wind particles which are subsequently analyzed with a 
mass spectrometer. The SWCE was one of the experiments flown 
on the Apollo 11 mission as well as several others and a sum¬ 
mary of the results from these missions is shown in Table 2-4. 

2.12 SWS-Solar Wind Spectrometer [Synder et al., 1970] 

The solar wind spectrometer was designed to investigate 
properties of the solar wind and its interaction with the moon 
It consists of seven symmetrically arranged Faraday cups, one 
directed vertically, and the other six at 60° angles from the 
vertical. This configuration allows observation of any solar 
wind anisotropies that may exist. Energy spectra may be 
obtained by modulating the voltage on the grids that precede 
the Faraday cups and noting the cup current changes. The 
SWS cycles through 14 energy steps ranging from 10 eV to 
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10.4 keV for protons and from 10 eV to 1480 eV for electrons 
in the high energy range and each energy step is reduced by 
a factor of 1.68 in the low energy range. 

Solar wind proton spectra obtained in this way show pro¬ 
ton velocities of typically 400 to 600 km/sec. The SWS does 
not generally see solar wind particles during the lunar 
night, indicating a plasma void behind the moon. 

The power for all of the experiments described, with the 
exception of the LPM and the SWCE is provided by a Radioisotope 
Thermoelectric Generator (RTG) which provides a continuous 
70 W. The individual experiments are connected to the ALSEP 
Central Station by Mylar coated ribbon-like cables where the 
data from each of them are collected and transmitted to earth 
through a high gain antenna mounted atop the central station 
and aimed at the earth. These data are received at various 
earth based tracking stations,recorded on magnetic tape and 
distributed by NASA to the various investigators. 
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Chapter 3 

This chapter is a description of the model used to study 
the behavior of suprathermal protons that are ejected into 
the solar wind in an upstream direction from the earth's bow 
shock. 

3.1 Equations of motion 

Consider a particle being ejected from the earth's bow 
shock with some initial velocity V p and define a coordinate 
system in which the z-direction is parallel to the interplan¬ 
etary magnetic field B, y is in the direction of -V sw X B 
and x completes the right-handed system. If the origin of 
this coordinate system is fixed at some point on the surface 
of the shock, then the particle sees an electric field 

"e = -V X*B, due to the convection of the interplanetary 
sw _i 

field lines past the earth at the solar wind velocity V sw . 

The equation of motion that describes the particle's 

trajectory upon leaving the shock is, 

m“P = q ( E + V X B } 3-1 

dt P 

in a coordinate system with B = Bz, one can choose E = Ey 
without loss of generality since E is always perpendicular to 
"b. The solution of equation 3-1 gives the time dependent 
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coordinates of the particle trajectory. 


v -V i 

x(t) _ — x °^ -JsinCOt + (1-cosCdt) + v fl t 


y (t) =—X2 sinCdt + 
CO 


(cosO)t-l) 


z (t) = v t 


3-2 


3-3 

3-4 


Equations 3-2 and 3-3 are the parametric equations of a 
cycloid in the x-y plane. For a derivation and discussion 
of these equations, see the appendices. The motion in the 
z-direction is uncoupled to motion in the x-y plane and hence 
if there is no initial velocity in the z-direction, the part¬ 
icle will move entirely in the x-y plane. 

If we transform the above equations to a frame of reference 
that is moving parallel to the x-direction with velocity V d , 
we see that the particle describes circular motion with center. 


mv 

__L° 

qB 




and radius 


(m/qB) 


xo d 

in the moving reference frame. That is 


. 2 2 

v.) + v 

A yo 


1/2 

3-6 

the guiding center 


of the circle moves in the x-direction with velocity V d =E/B, 
the drift velocity, (see figure 3-1) 
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Figure 3-la. Illustration of particle drift motion 
for a proton* in crossed E and B fields when the initial 
velocity of the particle is zero. 

Figure 3-lb. View of drift motion in the Y se -Z se P lane 
showing the effect of adding a finite z-component of 
the velocity. 

Figure 3-lc. View in the x se " Y se P lane showing the 
particle's resultant guiding center trajectory. 


/ 









Figure 3-1 



52 


Due to the uncoupling of the z-direction motion from the 
x-y planar motion, a simple expression can be written that 
determines the resultant motion of the guiding center of a 
that leaves the earth's bow shock with 
some intial velocity V p . The resultant velocity V r Q f the 
particle's guiding center is given by. 


3-7 


where V i s the component of the particle's initial velocity 
that is parallel to the interplanetary field and V d i s given 
by. 


V, —{EX B) /B 2 
a 


3-8 


3•2 Guiding center behavior 

Figure 3-2 is an idealized representation of the magnet- 
ospheric cavity as it exists in the free flowing solar wind. 

Xn a solar ecliptic coordinate system the interplanetary mag¬ 
netic field has a solar ecliptic longitude of approximately 
315°, the "garden hose angle", and lies more of less in the 
ecliptic plane. 

To understand the basic longitudinal behavior of the model, 



Figure 3-2. Diagram demonstrating the asymmetry of 
the E X B drift model and the basic guiding center 
behavior of a zero pitch angle particle. 
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consider a particle of zero pitch angle leaving the shock at 

point P in figure 3-2 and traveling along the field line with 

a velocity V , determined by its initial energy. If the in- 
11 

terplanetary field were stationary, the particle would move 
along the path determined by and it would intersect the 

lunar orbit at point B. Due to the convection of the field 
lines past the magnetospheric cavity at the solar wind veloc- 
ty, there exists an interplanetary electric field E= B 
such that the particle's guiding center follows a path that is 
determined by, 


(V,. 


B) B 


3-9 


and the particle will actually follow a guiding center path 
that intersects the lunar orbit at point A. If the moon is 
also at point A and the detector is looking in the right dir¬ 
ection, then the particle can be seen there if it is of the 
proper energy to be detected by the SIDE. 

Clearly, if a particle leaves point P with less energy or 
a higher pitch angle, i.e, a "shorter"V , the same drift vel- 

ocity vector added to its V will yield a resultant velocity 
ll 

vector that defines a guiding center path that intersects the 
lunar orbit at a point that is closer to the shock than point 
A. Conversely, a particle leaving the shock with greater par¬ 
allel energy, i.e. a "longer" , will have a resultant 


1 


I 


I 


f 
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velocity vector that will define a guiding center path that 
intersects the lunar orbit at a point further away from the 
shock than point A. 

That is, whether or not a particle of a given energy 
is detected at the lunar orbit is dependent upon the orienta¬ 
tion of the interplanetary B field, the magnitude and direc¬ 
tion of the solar wind and the particle's initial energy and 
pitch angle upon leaving the shock. 

From the above description of the energy selection mech¬ 
anism that is operating in the post shock region of the 
lunar orbit, it becomes clear that for a particle leaving 
the shock with a given parallel energy E , and hence a given 

_i 

initial parallel velocity V p , that there exists a point Q on 

n 

the bow shock that defines an energy cutoff point in the lunar 
orbit. If the moon has advanced in its orbit beyond the point 
C in figure 3-2, then particles of energy E-^ leaving the 
bow shock can no longer be seen by the detector. 

The cutoff point varies with energy and interplanetary 
field conditions, but one expects that in the oost shock cros¬ 
sing region of the lunar orbit, that the lower energy particles 
will begin to disappear from the data as the moon advances 
toward instrument nightime. This effect is observed in the 
data and will be discussed again in the chapter concerning 
the data, (chapter 5) 
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A completely analogous mechanism is operating on the inbound 
side of the bow shock and the appropriate vector diagram for 
that region is also illustrated in figure 3-2. Due to the 
geometric asymmetry in the field configuration, the data from 
the inbound side of the shock will show a somewhat different 
effect. The higher energy particles will be seen further away 
from the shock than will the lower energy ones, but none will 
be seen as far away from the shock as on the outbound side for 
a given field configuration. This is in agreement with the 
"dusk-dawn" asymmetry observed by Asbridge et al. [1968]. 

The above discussion required that the solar ecliptic lat¬ 
itude of the field be approximately zero since the moon's 
orbit lies approximately, (+ or - 5°), in the ecliptic plane. 
This would suggest that if the interplanetary B field were 
not in the ecliptic plane or close to it, then particles leav¬ 
ing the bow shock might not be seen by the detector. Indeed, 
if the field is far away from the ecliptic plane, particles 
leaving the shock along these lines could likely miss the 
lunar orbit entirely, (see figure 3-3) 

In the above discussion, the interplanetary field was 
taken to be approximately in the ecliptic plane. If the 
solar wind vector also lies approximately in the ecliptic 
plane, this requires that E=- v s ^. B be perpendicular to the 
ecliptic plane and hence will also lie in the ecliptic 
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Figure 3-3. Illustration in three planes of the dep¬ 
endence of the resultant guiding center trajectory on 
the interplanetary field latitude. 





But if B 


plane as will the resultant velocity vector V r . 
is not in the ecliptic, the resultant velocity vector will 
be out of the ecliptic plane. As can be seen in figure 3-3, 
whether or not a particles's guiding center trajectory inter¬ 
sects the lunar orbit is fairly sensitive to the interplanet¬ 
ary field latitude. An analysis of the data shows this effect 
and will be discussed in chapter 5. 

It should also be pointed out that there is an interesting 
dusk-dawn asymmetry in the direction of the interplanetary 
E field when the B field is at high positive or high negative 
latitudes. Figure 3-4 shows that for a B field latitude of 
+90° on the dawn side, that the interplanetary electric field 
points outwardly from the shock and similarly a B field lat¬ 
itude of -90° on the dusk side generates an E field that also 
points away from the shock. This gives rise to another mech¬ 
anism, as discussed by Manka [1972], for getting particles 
into the detector that has a definite dusk-dawn asymmetry 
associated with it. 

The foregoing description of the model has ignored any 
influence of the perpendicular energy of the particle and 
hence the size of its gyroradius. To gain understanding of 
this influence, the actual particle trajectories will now 
be discussed. 
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Figure 3-4. Diagram illustrating the dusk-dawn asym¬ 
metric interplanetary electric field E that occurs at 
high positive and negative values of the interplanetary 
B field latitude. 
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3.3 Actual particle tranectord.es 


The gyroradius of a 3 keV proton in a 5T field is about 
1500 kin which is comparable to the lunar radius (^ 1700 km). 
This implies that a strict interpretation of the above treat¬ 
ment, where the particle follows the guiding center path det¬ 
ermined by v r , is not quite correct. The actual trajectories 
depend on the particle’s initial velocity and pitch angle a 
and can be calculated from the parametric equations 3-2,3-3, 
and 3-4. The pitch angle is the angle between the initial 
velocity vector and the direction of the field which implies 
that 


3-10 


where 


and 


V = /v 2 + v 2 + t 2 ) 1/2 = (2 £/ m ) 1/2 3-11 

o 1 xo yo 70 / 1 ^ ' 

V - V sin a = ( V 2 + V 2 ) 1//2 3-12 

v o x 11 w y v xo yo j 


For typical conditions, V sv p400 km/sec, B=57 , several 
cases of proton motion are shown, to scale, in figures 
3-5 through 3-9. It is evident from these figures that part¬ 
icles whose guiding center trajectories intersect the moon 
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Figure 3-5. Scaled illustration of cycloidal motion 
of typical bow shock protons. Top curve is the case 
of all parallel energy and the bottom curves are cases 
of 1 keV of perpendicular energy distributed in three 
different ways. and are the values of the part¬ 

icle's initial velocity in the x and y directions res¬ 
pectively in units of km/sec. 

Figures 3-6 through 3-9. Additional cases of cycloidal 
motion for perpendicular energies higher than 1 keV. 

Note the interesting case in figure 3-6 where the x-comp- 
onent of the initial velocity is approximately equal to 
the drift velocity (283 km/sec). (See equation 3-6) 










£ =1 keV o( =•£ C° 



Figure 3-6 



£=2 kev 


°( =90' 



Figure 






Figure 3-8 




Figur° 3-9 



can enter the detector from a wide range of angles. The 
implications of this on the interpretation of the data are 
discussed in chapter 5. 
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3.4 Asbridge Model 

The model proposed by Asbridge et al. [1968] is very 
similar to the present treatment and the basic idea behind 
their model is illustrated in figure 1-6. The resultant 
velocity vector of the particle is determined by, 


where, is the particle’s initial velocity parallel to 
the field. This differs from the expression for the re¬ 
sultant velocity vector given in equation 3-7 by the factor, 

- ( V sw . B ) ^ 3-14 

which is just the component of the solar wind velocity par¬ 
allel to the field. When V sw is'perpendicular to B, the two 
expressions are equivalent. 

The basic difficulty with eqation 3-13 is that it appears 
to give the resultant guiding center velocity vector in the 
frame of the bow shock for a particle that is injected into 
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the solar wind along the B field in the shock frame. There 
is no question that if a particle left the bow shock in the 
direction defined by V f that it would continue along this 
path undeviated by any electric fields since 

■^ = ( E x + V f X B ) = 0 3-15 

nonrelativistically. However, if the particle actually had 
a velocity v j_ in the shock frame as indicated in figure 1-6/ 
it would not go in the direction defined by V f as indicated 
by the Asbridge model. 

This point can be clarified by considering the transform¬ 
ation between a frame of reference that is moving with the 
solar wind and the rest frame of the bow shock. Since V gw «c, 
the transformation is Galilean, 


y = y 1 3-16 

z = z' 

where the moving frame is being translated along the x-axis at 
the solar wind velocity. This implies that. 
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*V =~V' + V 3-17 

v sw 

If V' i s directed along *B in the moving frame, then the 
velocity seen in the rest frame is V. Equation 3-13 is actual 
ly just a transformation from a frame moving with the solar 
wind to the rest frame of the shock, i.e., 

V i = V* 3-18 

and 

V f = V 3-19 

It is not realistic to assume that one knows the inital 
velocity of the particle in the moving frame and therefore 
equation 3-7, derived in this chapter, is the appropriate 
expression to determine the resultant guiding center trajec¬ 
tory of a particle in the rest frame of the shock. It is 
interesting to note that whenV sw is not perpendicular to B, 
equation 3-9 yields a resultant guiding center trajectory 
that allows particles to travel to points at the lunar orbit 
that are more sunward than those predicted by equation 3-13. 

3.5 Conclusions 

The discussions in this chapter lead to the following 


conclusions. 
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a. ) There are particlular interplanetary field longitudes 
that are favorable to the observation of particles of a given 
energy at a given point in the lunar orbit and this fact gives 
rise to energy cut-offs in the post shock region of the lunar 
orbit and energy cut-ons in the pre shock region. 

b. ) There are particular interplanetary field latitudes 
that are favorable to the observation of particles in the 
post and pre shock regions of the lunar orbit, viz. those in 
the neighborhood of zero degrees. 

c. ) There is an inherent dusk-dawn asymmetry in the model 
that should be observable in the data. 

d. ) There is a dusk-dawn asymmetric interplanetary elec¬ 
tric field mechanism for getting particles into the detector 
that is dependent on high positive and high negative values 
of the interplanetary B field latitude. 

e. ) Although the Asbridge model of guiding center tra¬ 
jectories gives the correct qualitative picture of particles 
leaving the bow shock and the correct quantitative picture if 
interpreted correctly, the E X B model proposed in this thesis 
is more appropriate for the reasons discussed in this chapter. 
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Chapter 4 

-In light of the foregoing discussion, a computer model 
for the e: X "b mechanism has been generated. The description 
of this computer model and its connection to the physical 
model of the previous chapter will be given below. In order 
to generate a qualitative picture of what the detector should 
see as a function of time, it is necessary to calculate the 
components of a vector ^(t) that represents the look direction 
of the detector as a function of time and use this vector in 
conjunction with the vectors V ^(£,a ), the resultant guiding 
center trajectory vector/ and r, the position vector of the 
moon to determine whether or not a particle is seen. By de¬ 
fining a source with certain characteristics and "acceptance 
criteria" for the detector, a spectrum of counts vs. energy 
can be obtained and compared to the real data. What follows 
is a description of the computer model. 

4.1 Orbit 

In order to determine whether or not a particle leaving 
the bow shock will hit the moon, it is necessary to determine 
the position vector of the moon, r. Figure 4-1 shows the or¬ 
ientation of the lunar orbit with respect to an inertial sys¬ 
tem whose x-y plane is in the plane of the ecliptic. The or¬ 
bital inclination, Q , is exaggerated in the figure for clarity 
it has a value of about 5°. The x’ axis intersects the lunar 
orbit at perigee (denoted by 7r ) in the figure. 
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Figure 4-1. Illustration of the relation of the lunar 
orbit to the solar ecliptic plane showing time depend¬ 
ent angles relevant to the calculation of the position 
of the moon, n denotes perigee. Q is exaggerated in 
the figure for clarity. In reality q 5°. 
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and precesses with a period of approximately 8.5 years in a 
counterclockwise direction when viewed from a point on the 
z'-axis above the x'-y* plane. The line that divides the time 
between when the moon is above the ecliptic plane and the 
time when it is below the ecliptic plane is known as the line 
of nodes and it precesses in a clockwise direction when look¬ 
ing in the negative z-direction from a point above the eclip¬ 
tic plane. The line of nodes precesses with a period of ap¬ 
proximately 18.6 years. 

The moon moves in a slightly elliptical ( €~.05) orbit 
that lies in the x'-y’ plane with perigee occuring on the 
x’ axis. The distance from the earth to the moon as a func¬ 
tion of time is given by 


a (1- e 2 y __ 

1 + e cos ( <d t- t/» ) 


where, a = semi-major axis, € = eccentricity, 0)= angular 
velocity of the moon based on the sidereal period, $= argu¬ 
ment of perigee, and t begins at ascending node. 

The angles, relevant to the Eulerian transformation, as 
functions of time are determined by, 


0(t) = (p - ft l on t implies CW precession 

0 o = argument of line of nodes at t=0 4-2 

=18.6 yrs. 


0 = 2 i r / T 

Ion 


Ion 


Ion 
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«|» 0 = argument of perigee at t=0 
= 8.6 yrs. 


(1, = 27t/ t. 

loa loa 


6 (t) — const. 


loa 
6= 5° 


The transformation from a coordinate system, S’ 1 , that is 
fixed in the moon, to the solar ecliptic system, S ge is accom¬ 
plished by first transforming from S'’ to S* and then from 
S' to S and finally from S to S se . That is if _A, represents 
the matrix transformation from S' 1 to S', J3. represents the 
transformation from S' to S and _S_ represents the transforma¬ 
tion from S to S se then, 

~r' = jk~r ' ' 

r=_Br' 4-4 

"r = _S.Tr 
se — 

which implies that, 

~r = S B _A.Tr' ' 
se -- 

The coordinate system S t! is simply rotating with respect 
to S' so that the angle between the x''-axis and the x’-axis 
is, 

0= Cii t— t/f(t) 4-5 
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cos© -sin© 0 \ 

sin© cos© 0 4-6 

0 0 1 J 

The transformation from S' to S is given by equation 4-7, 
shown on the next page, 

The inertial system S is stationary with .respect to the 
fixed stars but the x-y axes of the SE system rotate slowly 
with respect to the inertial system due to motion of the earth 
in its orbit about the sun. That is, the x ge ~axis rotates with 
respect to the x-axis about a degree per day in the counter¬ 
clockwise direction as viewed from above the ecliptic plane. 
Therefore, the transformation from the inertial system to the 
solar ecliptic system is given by the matrix, 


and since z xl f z' r 

A = 


where, 


f cos £ -sin £ 0 

S = | sin £ cos £ 0 

0 0 1 


£ = £ + (2 TT t) / T e 


t = sidereal period of the earth 
e 

Therefore the components of the vector r in the SE system 
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where r 



C S B A r” 
( S B A )y X r x 
< S B A } zx r x' 


and the matrix elements of SBA are shown on the next page. 


4.2 Detector hook Direction 

If one know the orientation of the detector look vector, 
"qi i » a coordinate system, S' ' ’ , that is fixed at the ALSEP 
site, then the components of this vector can be expressed in 
a S.E. system, S^, by making an Eulerian transformation 
from S’'’ to S ge . 

The moon's rotational axis makes an angle with the 
ecliptic plane of about 1.5 degrees and the rotational period 
is equal to the sidereal period. Figure 4-2 shows the rela¬ 
tionship between the coordinate system fixed at the ALSEP 
site and the solar ecliptic system. 

The Apollo 12 and 14 instruments were deployed very near 
the lunar equator and for the purposes of this treatment we 



(SBA) 


xx 


(SBA) 

yx 


cos £ j”cos0 (CGS^COS^. - cos0 sin<£sintfr)-sin© (sintfr costf. + cos0sin<£cos$] 
-sinf [cos 0 (cosV'sin£ + cos0 cos<£ sin iff) -sin 0 {sinifr sin<£ - cos 6 cos0 cos iff )j 

sin^ [cos© {cos iff cost}} - cos Q sin 0 sin iff)- sin© (sin^cos0 + cos 6 sin<£cos^)J 
+cosf jcos © (cosiA sim£ + cos 9 cos 0 sin -sin 0 (sin« sin<#>-cosfl cos^cos^.)] 


(SBA) 


zx 


(cos© sin^ + sin© cos iff ) sin# 


Equations 4-9 


-j 



Figure 4-2. Relationship of the ALSEP coordinate sys¬ 
tem S' 11 , or system 1, to the SE coordinate system. For 
deployment at the sub-earth point, x ,M points to earth 
and remains approximately that way due to the equality 
of the lunar sidereal period and the rotational period, . 
Small changes are due to physical librations. 
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will assume that they were deployed exactly on the equator. 
Then at the ALSEP site the unit vector x' 11 is in the direc¬ 
tion of the local vertical, % ' 11 points north and in the case 
of equatorial deployment, z ,,f is parallel to the rotational 
axis of the moon. §'' 1 completes the right-handed system 
and is directed eastward. S rr% - system 1. 

The instrument, l assume Apollo 14), was deployed at the 
ALSEP .site so that the body axes were coincident with x'' 1 , 

■y’ ’ 1 and 1 ', The vector that represents the.look direction 
of the detector in the S’'’ system is then, 

D cos 15° 

D"’ = Id sin 15° 

\ 0 

The Eulerian transformation from S' 11 to S ge is given on the 
following page. Using this transformation D in the SE 
system is, 



that is 


*D = Td' 

' 1 

D x ^ ^ XX D x 

' + ( T 

D y = <t ) yx c 

' +. ( T 

D Z - < T 

' + ( T 


4-12 


4-13 






= cos£ (cost? cosA - cosy sin \ sin t?) -sin^ (cost? sin A. + cos y cos A sinr? ) 
T ~ sin^ (cost 7 cosA - cos y sin A sinr? ) +cos^ (cost? sin A. + cos y cosA sin 17 ) 


T - siny sin T? 


= -cos^ (sinT? cos A + cos y sin A cost? )-sin£ (cost? cos A cosy -sin A sinr? ) 
— -sin^ (sini? cosA + cosy sin A cost? )+cos^ (cost? cosA cosy -sin A sinT? ) 
T = siny cost? 

yz 

T = cose siny sin A + sine siny cosA 
zx 3 3 

T zy = sin£ siny sinA - siny cosA gos^ 

T :— cos y 


Equations 4-11 


At this point, it is possible to calculate the resul¬ 
tant velocity of particles leaving the bow shock, V r , the pos¬ 
ition vector of the room r, and the detector look direction 
all in a solar ecliptic coordinate system. 

4.3 Shape of the shock surface ; 

The source of the particles in this computer model is an 
idealized mathematical representation of the magnetospheric 
cavity that emits particles with a flat energy spectrum, 
isotropically at each point on its surface. The mathematical 
shape is a paraboloid of revolution that is fitted to the 
true shape of the bow shock surface. This fitting procedure 
will be discussed below. The model retains all of the 
qualitative features of the real bow shock and at least two 
of the quantitative features, the distance of the nose to the 
center of the earth and the distance from the center of the 
earth to the intersection of the lunar orbit with the shock. 
The paraboloid is shown in figure 4-3 and the equation in 
solar ecliptic coordinates is, 


x = -d ( y 2 + z 2 ) + d 4-14 



where d is the distance from the nose of the shock to the cen¬ 
ter of the earth and D is the distance from the center of the 
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Figure 4-3. Idealized bow shock surface represented 
by the equation X 2 =-d/D 2 (Y 2 +Z 2 )+d in a SE coordinate 
system. 




«* 



Figure 4-3 
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earth to the shock, measured along the y-axis, The shock 
stand-off distance d, can be calculated by equating the solar 
wind kinetic pressure to the magnetic pressure provided by 
the earth’s dipole field so that 


and the parameter D, can be determined from knowledge of the 
positions of actual bow shock crossings, cf. fig. 4-4a. 
That is by knowing the time of the bow shock crossing, the 
S.E. longitude can be obtained from the ephemeris data. The 
time of crossing can be determined accurately from the SIDE 
data. An example of this is shown in figure 4-4b. Solving 
equation 4-14 for D and neglecting any z dependence. 




D was calculated in this manner, but a constant value for d, 
(/s/12 P. ) was chosen over using the expression 4-15. 
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Figure 4-4a. Diagram illustrating the technique for 
fitting the idealized shock surface to the observed one. 
Knowing the distance d and measuring and Y^g allows 

determination of D from equation 4-14. 

Figure 4-4b. Illustration of how the postion of the 
bow shock is accurately determined for use in calculat¬ 
ing D as discussed above. 
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Figure 4-4 






4.4 Particle emission 

Each point of the surface of the shock is taken to be an 
isotropic radiator of particles with a flat energy spectrum. 
That is, the same number of particles are emitted for each 
energy and the spectrum looks the same in all directions from 
each point of the surface. The total number of particles 
emitted at a given pitch angle a is proportional to the area 
on a unit sphere between - ^ and-j-^ . See figure 4-5a. 

The number of particles available for detection is there* 
fore weighted toward higher pitch angles by a factor of sin a 
If n is the number of particles per square centimeter passing 
through a unit sphere centered at the point, then the number, 
dN, emitted at a given pitch angle a is. 


dN = ndA 

dN = 2jt n sina da 
dN = const, sina 

for a constant interval da . 

The total number N, emitted over the hemisphere, center- 
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Figure 4-5a. Shows that for an isotropic radiator 
the number of particles emitted at pitch angle ©< 
depends on . 

Figure 4-5b. Illustrates that all of the source points 
in the field of view of the detector are tested con- 
stituing a numerical integration over the source to 
determine the counts observed at the detector. 





83 


In tiie computer model , n was chosen to be of the same order 
of magnitude as observed in the data and pitch angle inter¬ 
vals, Act t of 10° were used. T’he initial values of the ener¬ 
gies ranged from 10 eV to 3500 eV. 

4.5 Acceptance Criteria 

£here are essentially two criteria to consider in deter¬ 
mining whether or not a particle is to be counted by the de¬ 
tector on the lunar surface. 

1. ) the particle must*come from the direction in which 
the detector is looking and 

2. ) the particle must have a source point on the sur¬ 
face of the shock, 

* (it is possible to relax this criterion under certain cir¬ 
cumstances; see discussion in chapter 51. Due to the finite, 
size of the detector look cone, there is more than one dir¬ 
ection from which the particles can come and still enter the 
detector. If 5 is the angle associated with the look cone, 
then since determines the direction of the guiding center 
path and “D determines the look direction of the detector, 
a particle can enter the detector if the following relation 
is satisfied. 


A A 


V ■ D ^ cos s 


4-19 
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The yector ( determines the direction in space of the 
guiding center trajectory and the vector*r determines the 
point in space through which the particle must pass. If the 
line, determined by V^. and the point *r, intersects the sur¬ 
face of the shock, then a source point exists. The equation 
of a line in three dimensions, passing through the point r, 
with direction V , is. 


4-20 


The simultaneous solution of this equation and the equation 
of the paraboloid of revolution eq. 4-14, yields a quadratic 
in any of the variables x,y, or z. Then if the discriminant. 


C B 


4 A C ) 


1/2 


is real, that is if 


a 4 4 A C 


a source point exists. In the case of a quadratic in y, 


^ V r^ d V 2 {V rs /V ry ) r V rz /V ry ) r y 


+■ (v /V ) 2 

\ rz ry/ 


)[(V /V } r - r 1 
I rz ry y z J 


2 2 2 

c * D (r /d) + (V /V ) r - 2 r r (v /v ) 
x rz' ry y y z rz' ry 

-D 2 - (D 2 r V /d V ) 
y rx' ry 
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It is now possible to outline the scheme for generating 
data with the computer model. 

a. ) the moon is advanced in the orbit by incrementing 
the time, thus determining the position“r(t), and the detec- 
tor look direction D(t), 

b. ) at each position, the values of a ) are cal¬ 

culated based on B and . 

sw ^ 

c. ) the tests V r ,cz)*?^cos8 , and B > 4 A C are per¬ 
formed and if they are both satisfied, dN (£,<z ) is stored 
in the appropriate energy channel. 

To allow for the higher band pass at higher energies, 
all particles from energyto (*¥ that satisfy the 
conditions described above, should be counted in the energy 
channel £ . Due to the large amount of computer time requir¬ 
ed to test a nearly continuous range of energies, only the 
twenty center energies were used and each dN (£, a ) was 
multiplied by 10% of the center energy. This is in keeping 
with the qualitative nature of the model. 

Hidden in this technique is a numerical integration 
of the source points over the surface of the bow shock. That 
is, by making the above tests at each position r, over a wide 
range of energies and pitch angles, each source point that 
contributes a particle that satisfies the tests is counted. 
This can be seen by examining the special case shown in figure 


4-5b. 
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Chapter 5 

5.1 Introduction 

Two basic approaches were followed in the analysis of 
particle data and interplanetary magnetic field data to deter¬ 
mine whether or not the physical processes, discussed in 
chapter 3 are in operation in the pre and post bow shock re¬ 
gions of the lunar orbit. 

The first approach consisted of searching for the ex¬ 
pected correlation in interplanetary field configuration and 
observability of particles with the detector on the lunar 
surface while the second approach consisted of attempting 
to generate synthetic SIDE data, using the computer model 
and to discover in which ways these data resembled the real 
data. 

Simultaneous particle and magnetic field data were 
available from the Apollo 14 SIDE and Explorer 35 magnetome¬ 
ter respectively, for three periods of interest during May 
and June of 1972. The magnetometer data tapes were provided 
by D.S. Colburn, Ames Research Center, Moffet Field, Califor¬ 
nia, The regions in the lunar orbit corresponding to these 
perio,ds and their associated times are shown in figure 5-1. 
Discussions of the analysis and results of the analysis of 
these data follow under the subtitles of Real Data and 
Synthetic Data. 
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Figure 5-1. The three regions of the lunar orbit 
in which magnetometer and ion data were collected 
for the analysis done in this thesis. The distance 
between the tick marks represents the distance 
traveled by the moon in 1 day. 


Lunation 17/ 1972 




REGION from 

INBOUND 141 d 16 h 2l m 

_m 

TAIL 148 0 1 

OUTBOUND 153 d 19 h 24 m 


to 

143 d 16 h 2i m 

150 d 0 h l m 

156 d 19 h 24 m 


141=May 20 
148=May 27 
153=June 1 


Figure 5-1 
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5.2 Real Data 

In order to determine whether or not any correlation 
between B field configuration and particles observed exists 
in the data, it is necessary to devise a technique of com¬ 
paring the known B field components with the known particle 
accumulations of the SIDE for the same time interval. 

The obvious thing that comes to mind is to simply plot 
the component of the field of interest and the integral par¬ 
ticle flux for the same time interval and visually inspect 
these plots for any correlation. An attempt at this proce¬ 
dure was made and the result is shown in figure 5-2. Al¬ 
though a weak visual correlation is evident from these data, 
it is obvious that no concrete conclusions can be drawn 
from this type of analysis. In order to overcome this problem, 
we could divide this plot into small time intervals and con¬ 
sider each interval separately, noting the value of the field 
latitude and the number of particles observed at that lati¬ 
tude for each interval. That is, we would be interested in 
constructing from the data, a plot of the average number of 
particles observed when the interplanetary B field has a 
particular configuration vs. the field configuration. 

Such an analysis was done using the SDS 910 computer 
and data tapes from the Apollo 14 SIDE and the Explorer 35 
magnetometer. 
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Correlations in solar ecliptic latitude of the field 
and solar ecliptic longitude of the field, separately, were 
sought as were correlations with the field latitude and 
longitude considered jointly. Details of the computer analy¬ 
sis can be found in the appendix on computer programs. 

5.3 Results from the Real Data 

The results of the analysis using the real data are 
shown in figure 3 through 12 of this chapter, for the three 
periods of interest, which will henceforth be referred to as 
the inbound, tail and outbound regions. The inbound region 
corresponds to the pre-shock crossing region and the out¬ 
bound region corresponds to the post-shock crossing region 
of the lunar orbit. The inbound region consists of two days 
of data taken prior to entering the magnetosheath, beginning 
at 141 d 16 h . The tail region begins at 148 d 0 h and continues 
for two days? this period is entirely in the geomagnetic tail 
and the analysis in this region was done as a check on the 
validity of this technique. The outbound regions begins a- 
bout 3 hours after shock crossing at 153 d 19 h 24 111 and ex¬ 
tends for three days toward the sunset terminator. 

The format of the first three graphs is designed to show 
the average number of counts observed at a given solar eclip¬ 
tic latitude of the field vs, the latitude and the frequency 
of occurence of that latitude in units of, number of times 



90 


Figure' 5-2. Visual correlation of field latitude 
and integral flux. A correlation is especially 
apparent in the second half of this event. 

Figure 5-3. Correlation of average integral counts 
and SE latitude of B over the inbound region. Bot¬ 
tom part of the graph shows the number of times a 
particular value of the latitude was encountered 
in the two day period over which this analysis was 
done. 

Figure 5-4. Lack of correlation in field latitude 
and counts in the tail region is apparent in this 
figure. 

Figure 5-5, Correlation of field latitude and counts. 
Note the broad central peak indicating that particles 
are preferentially observed when the field is near 
the ecliptic plane. The enrichment of the high 
positive latitudes of this distribution is due to 
lunar atmospheric ions. 
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observed per number of days of observation vs, the latitude. 
That is the lower portion of each graph displays simply the 
number of times a particular value of the field latitude was 
observed during the two or three day period of observation 
and the upper portion of each graph displays the average 
number of particles observed when at that latitude. Average 
integral counts will be taken to mean the average that is ob¬ 
tained by dividing the total number of particles accumulated 
at a particular latitude over the period of observation by 
the total number of times that that particular value of the 
latitude was observed during the same period. 

Due to the fact that the energy bandpass of the detector 
is proportional to the center energy, using counts instead 
of integral flux does not always provide a clear physical 
picture from the data. In view of this, the distributions 
of observed flux vs. field latitude were calculated for the 
three regions as a check on the validity of using counts for 
this analysis. The results are shown in figure 5-6. By 
comparing these results with those in figures 5-3 through 
5-5, it is clear that using counts does not substantially 
change the form of the distributions. Comparison of the two 
results for the inbound region does show a fairly large in¬ 
flation of the high negative latitudes in the integral flux 
representation that is not as apparent in the representation 
using counts. The reason for this is that the events that 
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caused the inflation have lower energies t £ < 500 eV) than 
those of the broad central peak. Using integral flux tends 
to amplify the lower energies due to the smaller bandpass 
at these energies. The particles found in the upper and 
lower ends of the distributions of figures 5-3 to 5-6 are due 
largely to lunar atmospheric ions that are accelerated by 
the interplanetary electric field and the effect of these 
particles on the present analysis is discussea more fully 
in chapter 6. 

The distributions shown in figures 5-7 to 5-9 are analo¬ 
gous to those in figures 5-3 to 5-5 where the solar ecliptic 
longitude of the field was the parameter of interest. It 
should be noted that the analysis was carried out here us¬ 
ing the uncorrected values of the field longitude. This 
has no bearing on the physical results, and in any physical 
interpretation of these results, the scale in longitude may 
be simply shifted by 28°, i.e. the corrected value of the 
longitude is obtained by adding 28° to the uncorrected value. 
For a discussion of the error in longitude and the technique 
used to correct for it, see chapter 2. 

The last three graphs of this series display the results 
of considering the latitude and longitude jointly. The eval¬ 
uation of the parameters in this treatment is over wider in¬ 
tervals, 10° latitude by 20° longitude, but the same proce¬ 
dure was used to determine the value in each two dimensional 
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Figure 5-6. Correlation in field latitude and int¬ 
egral flux. High latitude events are enhanced by the 
geometric factors, (see Appendix 3). 

Figure 5-7. Correlation of field longitude with counts . 
for the inbound region. The distribution exhibits a 
double peak but is not well defined due to the scarcity 
of data for this region. 

Figure 5-8. Unexpected correlation in field longitude 
and counts in the geomagnetic tail. Note the well de¬ 
fined peaks in the frequency' distribution and the differ¬ 
ence in height of the two major peaks in the counts dist¬ 
ribution . 

Figure 5-9. Correlation in field longitude and counts 
for the outbound region. The two major peaks are about 
180° apart. 

Figure 5-10 through 5-12. Correlations of field latitude 
and longitude considered simultaneously with average in¬ 
tegral counts. These plots contain basically the same 
information that is in figures 5-3 through 5-9. 
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interval. The average integral counts in each interval of 
the grid is represented by the number of lines in that inter¬ 
val where each line represents five or ten counts. The ab¬ 
solute value of any of the numbers is of no consequence and 
no attempt is made here to display the numbers in a rigorous¬ 
ly quantitative fashion. These graphs were also produced 
with uncorrected longitudes and the same correction angle 
may be applied to them as was done for the case involving 
only the longitude correlation. 

After generally surveying these results, the first im¬ 
portant conclusion that can be drawn is that there is indeed 
a correlation between interplanetary field parameters and 
particle observations by the SIDE at the lunar surface. This 
statement can be made on the assumption that if the number of 
particles observed by the detector did not depend on the 
configuration of the field, then the average number of parti¬ 
cles observed at any given value of a field parameter would 
not be very different from the average number of particles 
observed at any other value of the field parameter. In 
fact, that is the reason for carrying out the analysis in the 
tail region. The mechanism of the model is not expected to 
be operative in the tail region and a lack of correlation in 
field latitude can be seen in figure 5-4. Although the lati¬ 
tude of the geomagnetic tail field is sharply peaked about 
0° in frequency of occurence, the average number of particles 
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observed at any of the latitudes is not too different from 
that observed at any other latitude. A clear lack of corre¬ 
lation in geomagnetic tail field longitude was not observed 
and this somewhat unexpected result will be discussed in a 
separate subsection below. 

With the foregoing in mind, a correlation of interplane¬ 
tary field latitude is apparent in figure 5-5 for the out¬ 
bound region. A rather broadly peaked distribution is ob¬ 
served with the peak centered around -10 to -15 degrees, 
i.e. with the interplanetary B field near the ecliptic 
plane. Although the peak is not sharp, for reasons that will 
be discussed later, there is a tendency for more particles 
to be seen, on the average, when the field latitude has a 
value that is in the neighborhood of 0°, i.e. with 
-25° < B lat 4= +15°. Particles are still seen, however, at 
high positive values of field latitude and these are probably 
due to the outwardly directed interplanetary E field at high 
positive latitudes that accelerates the particles directly 
into the detector - . (see fig. 3-4) 

The results for the inbound region, figure 5-3, a^e not 
as well defined, due primarily to the shorter period of ob- 
servation> but again there is a broad peak that centers near 
0° latitude with some particles still being seen at high neg¬ 
ative values of the latitude. In the inbound region, if the 
field latitude has high negative values, the interplanetary 



96 


E field again points out into the detector on this side of the 
magnetosphere Csee figure 3-4) accounting for these high neg¬ 
ative latitude particles. 

At this point it is possible to conclude that the ob¬ 
servability of particles in the inbound and outbound regions 
of the lunar orbit is dependent on interplanetary field lat¬ 
itude and the dependence does roughly correspond to that ex¬ 
pected from the model. In addition, there is indication from 
the data that the asymmetric interplanetary E field at high 
positive and negative values of the B field latitude is pro¬ 
viding some of the particles observed by the detector in 
these regions. 

Figure 5-9 shows a correlation in field longitude for 
the outbound region with a fairly obvious doubly peaked dis¬ 
tribution, the first peak coming at about 84° uncorrected 
and the second smaller peak coming at about 264° uncorrected, 
a separation of 180° in longitude. One would also expect 
from the model a doubly peaked distribution with approximate¬ 
ly the same size peaks, separated by about 180° since the 
polarity of the field is unimportant in determining the re¬ 
sultant guiding center trajectory of a particle for any par¬ 
ticular value of the field latitude. For example, if the 
field latitude is say -10° while the longitude is 84° and 
+10° when the longitude is 264°, the same physical situation 
exists with respect to the model and the polarity is unim- 
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portant. If, however, the latitude is also -10° at the 264° 
position in longitude, a different physical situation doe 
exist and the guiding center trajectories will depend on the 
polarity of the field. Figure 5-13 illustrates this effect. 
Inspection of figure 5-12 reveals that the latitudes of the 
peaks in figure 5-9 are not too different. That is, differ¬ 
ent physical situations exist for the two peaks. The plane 
of the guiding center trajectories of case (1) in figure 
5-13 makes a smaller angle with the plane of the lunar orbit 
than does the plane of guiding center trajectories of case(3), 
possibly accounting for t* i favorability of observing parti¬ 
cles in cased), i.e. the higher peak at about 84° longitude. 

The data from the inbound region, figure 5-7, also show 
correlation in longitude of the field with the number of par¬ 
ticles observed, but again, the shorter period of observation 
results in a less well defined distribution. Despite this, 
the distribution also shows a double peak with a less pro¬ 
nounced difference in height. A less pronounced difference in 
latitude for the two peaks is found by examining figure 5-10. 
The two peaks occur at latitudes much nearer the ecliptic 
plane thus lessening the effect shown in figure 5-13. These 
differences are consistent with the ideas discussed above 
concerning the actual physical situation occuring for the 


two different peaks. 
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Figure 5-13. Views in two different planes of the phy¬ 
sically different situations ocurring for field config¬ 
urations in which the longitudes are 180° out of phase 
and the latitudes are equal. 





Cases 1 and 2 physically equivalent with respect to model 
Cases 1 and 3 not physically equivalent 

Cases 1 and 3 are favorable to observation of particles 
but case 1 is more favorable than case 3 according to the 
data in figure 5-8. 


Figure 5-13 
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The results of the longitude analysis show that there 
is correlation in observability of particles with the SIDE 
and the solar ecliptic longitude of the interplanetary B 
field and this correlation further supports the model dis¬ 
cussed in chapter 3. 

5,4 Interpretation of tail passage 

As mentioned above, figure 5-8, exhibits a doubly peak¬ 
ed distribution with the peak corresponding to the field of 
the geomagnetic tail pointing in an anti-solar direction at 
about a factor of 2 greater than the peak corresponding to 
the tail field pointing in the solar direction. This may not 
be significant, but since there is no obvious lack of correla¬ 
tion in longitude and number of particles observed, the phy¬ 
sical situation in this region should be examined more closely. 

At the beginning of the tail region, (14 8 d 0 h l m ), the 
solar ecliptic longitude of the geomagnetic tail field is 
fairly stable around 358°, corrected, or with the tail field 
approximately in the solar direction. At approximately 
148^ 5^, the field longitude changes abruptly and becomes 
stable about the value 178°, or the anti-solar direction and 
remains there throughout the two day period. The fact that 
no more abrupt changes in field longitude, i.e. neutral sheet 
crossings, are observed over this period can be understood by 
looking at a diagram of how the relative distances from the 
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ecliptic plane of the moon and the neutral sheet at the long¬ 
itude of the moon vary with time over the tail region at 
60 R g . The curve R in figure 5-14 represents the distance 
of the neutral plane from the ecliptic plane at the longitude 
of the moon, taking into account the seasonal offset and the 
diurnal variation? an off axis pivot point is also assumed 
due to the effect of the solar wind aberration angle of 4°. 
The curve d of figure 5-14 is simply the distance of the moon 
from the ecliptic plane at 60 R g as it moves toward higher 
longitudes. The sense of the geomagnetic tail field, in the 
region in figure 5-14 that is above the curve R, is in the 
solar direction and the field has the opposite sense in the 
region beneath the curve R. The motions in this diagram 
should be considered from the point of view of an observer 
on the earth, looking in the anti-solar direction. The moon 
spends a short time in the northern lobe region of the tail 
field and then crosses into the southern lobe for the rest 
of this period. 

All that can actually be concluded from the result in 
figure 5-8, is that there is a tendency for more particles to 
be seen, by the detector on the average, when the tail field 
points in the anti-solar direction, i.e. toward the detector, 
than when the field points in the solar direction, away from 
the detector. Note that the detector looks earthward in this 
region of the orbit. There is evidence here for a shadowing 



Figure 5-14. View from earth, looking in the anti-solar 
direction. Curve R represents the distance in R £ of the 
magnetic neutral plane from the ecliptic plane as a func¬ 
tion of time. Curve d represents the path of the moon 
with respect to the ecliptic plane for the same time per¬ 
iod. Note that the moon spends most of the two day period 
in the southern lobe of the earth's field. 





Figure 5-14 
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effect by the moon, yan Allen et al. I1S.7U , have suggest¬ 
ed a model by which thermal anisotropies in the solar wind 
can provide an asymmetric access of protons to the magneto¬ 
sphere. The basic idea of this model is that interplanetary 
field lines connect to geomagnetic field lines in the dis¬ 
tant magnetotail (~100 R e ). When the interplanetary field 
has a southward component and is directed away from the sun 
it connects to the tail field lines as shown in figure 5-15. 

It can be seen from examination of this figure that the pos¬ 
sible shadowing seen in the present analysis is consistent 
with the Van Allen model. By examining other features of 
the data in this region, for example the energy spectra, a 
better determination of the origin of this phenomenon could 
be made. 

5.5 Synthetic Data 

The basic philosophy underlying the second approach has 
been described in chapters 3 and 4 and the results of the ap¬ 
plication of that philosophy and a brief restatement of the 
scheme for generating the synthetic data are given in the fol¬ 
lowing sections. 

5.6 Brief restatement of the scheme 

a.) the position vector "r(t) and the detector look dir¬ 
ection D(t) are calculated. 





Figure 5-15. Van Allen model for asymmetric access of 
solar wind protons into the magnetosphere. The clustered 
arrows on the field lines represent proton thermal an¬ 
isotropies in the solar wind. [Van Allen et al. (1971)]. 
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b, l at each value of t, the resultant velocity vectors 
V r C £, a ) are calculated based on B and 

c. ) two constraints are then imposed to consider the par¬ 
ticle countable. 

1. ) the angle between V r ( £ , “1 an<a D * s required to 
be less than some limiting angle, the value of which is the 
subject of the next subsection. 

2. ) the particle must have a source point on the 

shock surface. 

5.7 Detector Look Cone 

The detector look cone will be taken to mean the solid 
angle centered on the detector look direction t and having 
half angle 6 . The nominal value of 6 is ~3° as deter¬ 

mined from the angular calibration of the instrument using 
a parallel ion beam in the lab. It must be kept in mind that 
when the detector is observing ions of large gyroradius com¬ 
pared to the dimensions of the entrance aperture of the in¬ 
strument, the look cone of 3° half angle does not prevent the 
detector from seeing particles whose guiding center traject¬ 
ories come from directions that are not in the field of view 
of the look cone. This is especially important in interpret¬ 
ing the data in this chapter since the gyroradii of bow 
shock protons are typically hundreds of kilometers and the 
guiding center trajectory does not represent the actual path 
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of the particle. One would expect that if the SIDE could see 
particles only in the 6° look cone that the distributions 
shown in figures 5-3 through 5-8 would have sharp peaks about 
6° wide. One would also expect that if the effective size 
of the look cone was say B= 90°, then when the detector was 
very near the shock it would always see particles since the 
angle subtended at the detector by the shock surface would be 
very large. The distributions are very broadly peaked in¬ 
dicating that the effective cone angle for observing bow 
shock protons is much larger than the nominal 6°, but that 
the half angle is not as large as 90°. 

The influence of this effectively larger cone angle on 
the computer model is to relax the criterion that the resul¬ 
tant velocity vector lie within a 3° half angle of the detec¬ 
tor look direction. It was determined from testing the com¬ 
puter model that a cone half angle of about 60° gave the best 
agreement with the real data. For further discussion of the 
detector look cone and its relation to a particle's pitch 
angle see the discussion of the range of observable pitch 
angles later in this chapter. 

5.8 Application of Computer model with idealized condi tions. 

The first application of the computer model was done 
using the Rice IBM 370 computer, using artificial input para¬ 
meters, with the exception of those parameters relating to the 
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lunar orbit. Explorer 35 data was not available at this 
time so typical values of the interplanetary field parameters 
and solar wind velocities were used to test the qualitative 
behavior of the model. The result, shown in figure 5-16, was 
generated using a solar wind velocity of 400 km/sec, a field 
latitude of 0° and a field longitude of 315°. Ninety pitch 
angles were tested at each energy and a spectrum was generated 
for every 40 minutes of lunar travel for the 5 day period 
shown. 

Particles with energies of less than about 500 eV were 
not observed and only the highest three or so energy channels 
(3 to 3.5 keV/q) contained counts when the last spectrum was 
generated. The cutoff shoulder moved steadily toward higher 
energies with increasing time for each spectrum some particles 
were seen. The behavior of the model exhibited in figure 5-16 
is roughly characteristic of what is actually observed dur¬ 
ing the post bow shock crossing regions one each lunation. 

The spectra from the real data just after bow shock crossing 
are generally peaked around 500 eV to 1250 eV and after 4 or 
5 days only higher energy particles are seen if any at all. 
Figure 5-17 shows a typical post bow shock region for compari¬ 
son. 

It should be noted that in this highly idealized case/ 
at least the gross features found in the real data are re¬ 
produced. The low energy peaks in the real data are due to 
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Figure 5-16. Synthetic data using typical interplanetary 
field conditions, namely B lat =0 a *id B long =315 °* T ‘* ie ^ 
energy channels of the instrument are plotted on the hor¬ 
izontal scale representing 10 eV to 3500 eV; log^g of the 
average counts is plotted vertically and time runs up the 
oblique axis. 
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Figure 5-16 




Figure 5-17. Typical SIDE counting rate spectra from the 
dawn-side bow shock region. 
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lunar atmospheric ions. 

The solar wind velocity and the interplanetary magnetic 
field are highly time varying quantities and only if this is 
taken into account in generating these spectra can one expect 
to effectively test the model of chapter 3 by generating 
synthetic data, A more elaborate procedure for generating 
synthetic data would be to calculate actual particle traject¬ 
ories which may be done in future analyses of this problem. 

5.9 Inclusion of interplanetary magnetic field data. 

When Explorer 35 magnetometer data became available for 
the outbound region, discussed above, the computer model 
was refitted for the SDS 910. Ephemeris data tapes for 
Apollo 14, containing information on lunar position and de¬ 
tector look direction also became available so that it was 
no longer necessary to calculate the orbital parameters as 
described in chapter 4. The field information from the mag¬ 
netometer and the position information from the ephemeris 
were synchronized and fed into the model and analyzed at 8 
minute intervals using the scheme discussed above. The solar 
wind velocity was still taken to be 400 km/sec throughout. 

The counting rate spectra from the computer model, at 
8 minute intervals, for the outbound region are shown in 
figures 5-18 and 5-19. The gaps in time are due to the non- 
continuous data from the Explorer 35 magnetometer. It is 
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Figure 5-18, Synthetic data generated using the real 
values of the magnetic field parameters and the real 
ephemeris data as input to the program (see Appendix 3). 
Data gaps are due to the gaps in the Explorer 35 tapes. 

Figure 5-19. Continuation in time of figure 5-18. 

Figure 5-20. Real ion counting rate spectra for the 
time corresponding to that in figure 5-18. The high 
sharp peaks are transmission noise. 

Figure 5-21. Continuation in time of figure 5-20. 
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immediately obvious that the inclusion of the magnetic field 
data causes more structured synthetic counting rate spectra, 
but again it can be seen from examining these results that 
the observed particle energies are in the higher regions of 
the spectrum. That is, particles with energies less than 
about 500 eV are generally not observed. It is also interest¬ 
ing to note that the spectra are no longer always continuous 
in energy after reaching the cut-on energy; multiply peaked 
spectra are readily apparent. The results shown in figures 
5-18 and 5-19 resemble real counting rate spectra much more 
closely than those of figure 5-16, even without including 
the changes in solar wind velocity. 

The real counting rate spectra for the outbound region 
are shown in figure 5-20 and 5-21 for comparison. Although 
there are some obvious differences, the data generated by the 
computer model are not completely unlike the real data with 
both exhibiting structured spectra with counts mostly in 
the higher energy channels. The noisy spectra observed in 
the early post bow shock region of this lunation is not al¬ 
ways present and could be due to some turbulent effect, seen 
only while the moon is close to the shock itself or possible 
temperature effects on the instrument are responsible for 
this noise. Unfortunately, event by event comparisons can¬ 
not readily be made due to the time gaps. 
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After generating the data shown in figures 5-18 and 
5-19 r an analysis, using the real data, was applied and this 
result is shown in figure 5-22. Due to the fact that the 
spectra were generated only every 8 minutes, the number of 
times a particular value of the field latitude was observed 
is considerably less than the number observed in the analysis 
of the real data and hence, the distribution of average in¬ 
tegral counts is not as well defined. The results for the real 
data for the same period are shown in figure 5-23 for compar¬ 
ison. 

It may seem circular to analyze the generated data under 
the philosophy for the first approach since the model was 
forced to operate in generating these data whereas in analyz¬ 
ing the real data, we were looking for evidence that the mo¬ 
del was operating. It is, nevertheless, interesting to note 
that by forcing the model, a distribution is obtained that 
is very much like the actual distribution. 

5.10 Range of observable pitch angles. 

Due to the rather narrow field of view of the detector, 

{3° half angle), and the finite extent of the bow shock sur¬ 
face, some particles leaving the shock surface with large 
pitch angles may not be seen by the detector. It is possible 
to estimate the range of observable pitch angles at any given 
position in orbit by considering the simple geometry of fig- 



Figure 5-22. Correlation in field latitude and counts 
using synthetic ion data, as input. 


Figure 5-23. Correlation of field latitude and counts 
for the same time period using real ion data. 
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ure 5-24. The actual bow shock surface is approximated by a 
flat surface subtending approximately the same angle at the 
detector that the actual surface would subtend. The points 
L represent the limits of the source region, D represents 
the radius of a cross section of the bow shock in a dawn- 
dusk meridional cut and R is the geocentric distance of the 
moon. S is the half angle of the detector look cone, y is 
the angle subtended at the moon by the bow shock and a 's 
represent pitch angles. 

It can be seen from the figure that a particle leaving 
the shock at point C must have a pitch angle of S or less to 
be detected by the instrument on the lunar surface. However, 
if the particle comes from the point L it must have a pitch 
angle in the range 

y-S< a < y + s 

to be detected. Although each point on the surface defines 
such a range of pitch angles, the particles coming from the 
points L define the upper limit of the pitch angle range 
for a given distance from the shock. That is, particles 
with pitch angles 

a > y + S 

cannot get into the detector from any point on the shock at 
the distance defining y . 



Figure 5-24. Illustration of the relationship between 
the detector look cone angle, the angle subtended at the 
moon by the bow shock and the pitch angle of the particle. 
The figure shows the moon at last quarter in the dusk-dawn 
meridional plane. The dotted particle trajectory represents 
the upper limit to the pitch angle at this distance. 



to sun 



Figure 5-24 
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When the moon is very near the shock y is very large 
and the particles of high pitch angle may be observed while 
as the moon moves away from the shock y gets smaller, limit¬ 
ing the range of observable pitch angles. The limit of ob¬ 
servable pitch angles for a case in which the a-range is 
fairly restricted can be estimated by calculating Y when the 
moon is on the dusk-dawn meridian. For this case, 

D ~ 20 R 

e 

R /v 60 R 

e 

Y = tan -1 /20\ = 20° 

H 

a <23° 

This limit is probably one of the more restrictive due 
to the fact that as the moon moves further toward lunar 
nightime the detector sees a larger shock surface, thereby 
increasing y and as mentioned above, y is larger when the 
moon is nearer the shock surface. 

The effect of this limited range of c.tch angles was not 
considered in constructing the computer model, but the effect 
of this neglect should not be significant. Particles 
with high pitch angles have small II velocities causing their 
resultant guiding center trajectories to be more in an anti¬ 
solar direction than ones with small pitch angles indicating 
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that as one gets further away from the shock it becomes more 
difficult to see a high pitch angle particle based on consid¬ 
erations of guiding center trajectories alone. Rot reject¬ 
ing a particle in the model because of its high pitch angle 
does not significantly degrade the synthetic data since the 
particle would probably not be seen anyway because of its 
largely anti-solar guiding center trajectory. 

The above considerations neglect the fact that the 
pitch angle is not constant with regard to the guiding center 
trajectory, but similar behavior of the pitch angle limits 
can be expected if this effect is taken into consideration. 

5.11 Energy density of upstreaming protons 

To estimate the effect of the upstreaming protons on the 
incident solar wind it is appropriate to calculate the energy 
density of the particles and compare this with the energy 
density of the incident solar wind protons. 

By measuring the counterstreaming flux of the particles 
seen at the detector and assuming that each observed particle’s 
energy is due to its linear motion into the detector, one 
obtains an energy density of 

-13 3 

t? b s ~ 7 X 10 ergs/cm 

for the 1 to 3.5 keV particles seen on a typical dawn side 
bow shock region. The incident solar wind energy density is 
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of the order of 7 X 10“ 9 ergs/cm 3 which implies that 

\ s / ^ W ~ 10 ~ 4 (1 keV to 3 * 5 kev) 

One may conclude that the bow shock protons observed by 
SIDE would probably not cause a significant perturbation of 
the incident solar wind. 

Other investigators have calculated energy density for 
the counterstreaming protons that are considerably higher 
than this value, however none of these other values are 
greater than about 1% of the incident solar wind energy 
density for quiet conditions. Lin and co-workers [1974] 
calculate a value of 10 -12 ergs/cm 3 for protons in the range 
29 to 100 keV. Evidence from their data indicates that there 
is a sharp cutoff in counterstreaming protons above 100 keV 
which implies that 100 keV is probably the upper limit in 
energy for bow shock protons. Asbridge et al. [1968] calcu¬ 
late an energy density of 10 ^ ergs/cm 3 in the range of 3 
to 6 keV which is 1% of the incident solar wind energy 
density. Scarf et al. [1970] have observed values of the 
counterstreaming energy density that are comparable to the 
incident solar wind energy density at times. However, it 
seems that large scale, long time perturbations of the inci¬ 
dent solar wind by the counterstreaming protons is not a 
significant factor in determining the general behavior of 
these particles. 
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5.12 Conclusions. 

The results of the analyses presented in this chapter 
give rise to the follov7ing conclusions: 

a. ) the observability of particles at the lunar surface, 
during the pre and post shock regions of the lunar orbit, is 
at least partially dependent on the configuration of the in¬ 
terplanetary magnetic field. 

b. ) the evidence from the analysis of the data and 
from the behavior of the model suggests that the model dis¬ 
cussed in chapter 3 is probably operative in the pre and 
post shock regions of the lunar orbit and is responsible for 
the moderate to high energy particles observed in these re¬ 
gions. This result is not as well established in the pre 
shock region as in the post shock region and more data are 
needed to further analyze this region. 

c. ) The non-zero gyroradius of the particles cannot be 
ignored in the computer model but the effect it has on gen¬ 
erating synthetic data can be approximated by enlarging the 
detector look cone angle to include a larger range of guxd- 
ing center trajectories. 

d. ) the energy density of the counterstreaming protons 
is not sufficient to significantly perturb incident solar 


wind conditions. 


120 


Chapter 6 


6.1 Discussion 

There are possible sources other than the bow shock of 
moderate to higher energy protons, £ >750 eV, seen in the bow 
shock regions, viz. 

a. ) lunar atmostpheric ions 

1. ) interplanetary E field ions 

2. ) surface E field ions 

b. ) solar wind ions and suprathermal solar wind ions 

c. ) cosmic rays 

d. ) unknown sources such as lunar bow shock protons 
The "bow shock region" will refer to the regions of the 
lunar orbit extending from shock crossing to about a day 
or so past the sunset terminator on the dawn side and the 
analogous region before shock crossing on the dusk side of 
the magnetosphere. 

The reasons that it is unlikely that these sources can 
account for the observations completely will be pointed out 
below. 

a.) lunar atmospheric ions-As discussed by Lindeman [1973] 
and Manka [1972], ions from the lunar atmosphere can be ac¬ 
celerated into the detector by the- interplanetary electric 
field when the field lies in the ecliptic plane. By 
considerations of surface density and temperature, ion scale 
height, ionization rate, and strength of the electric field, 
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the expected energy spectrum observed by the SIDE is given by, 
F ( £ ) = ex P [ “ £/eEHJ 6-1 

where, 

n 0 =surface number density 
p=ionization rate 
e=charge on the ion 
H=scale height 



£ =ion energy 

Using appropriate values of these parameters for the 
Apollo 15 event observed on day 312, 1571 Lindeman [1973] 
has determined the predicted flux values and compared them 
to the observed fluxes. The result is shown in figure 6-1 
(after Lindeman [1973]). The flux drops off rapidly after 
about 250 eV indicating that the higher energy ions observed 
in this same region are not lunar atmospheric ions which 
have been accelerated by the interplanetary electric field. 
It should also be pointed out that the distribution obtain¬ 
ed in figure 5-5 indicates that most, of the particles were 
observed when the interplanetary electric field was not in 
the ecliptic plane in contrast to the results obtained by 
Lindeman [1973],(figure 6-2), and Manka [1972],(figure 6-3), 
that lunar atmospheric ions are seen when E. 


is in the 
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Figure 6-1. Predicted and observed energy spectra of 
lunar atmospheric ions. Note the sharp drop in energy 
after about 250 eV. [after Lindeman (1973)]. 

Figure 6-2, Integral flux vs. angle out of the ecliptic 
plane of the interplanetary*E field. The integral flux 
drops smoothly as the field goes out of the ecliptic. The 
bottom part of the graph displays £ /|f| vs. angle out 
of the ecliptic plane. [Lindeman, 1973] 

Figure 6-3, Correlation in integral flux and B^. 0 

implies an E field that is about in the ecliptic plane 
and pointing into the detector. [Manka (1972)] 
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ecliptic plane. 

Indications are that at the terminator the largest value of 
an attractive lunar surface potential is on the order of -100 V 
[Fenner et al., 1973 ], which implies that essentially thermal 
atmospheric ions could not attain energies much greater than 
about 100 eV before being detected by the SIDE. The lunar 
surface potential is obviously not responsible for the part¬ 
icles of energy greater than about 100 eV. It is difficult 
to imagine any other local electric fields of sufficient 
strength to accelerate atmospheric ions to the observed ener¬ 
gies . 

b.) solar wind ions-The thermal energy of a solar wind 
proton is on the order of 10 ev and its energy due to bulk 
flow motion, (V 400 km/sec), is > 1000 eV. The angle 
between the look direction of the detector and the solar wind 
flow in the bow shock region eventually becomes >9 0° # yet the 
high energy flows are still observed. 

Although local lunar surface magnetic fields as 
high as 300 gammas have been measured, Dyal et al. [1972], 
these fields are not of sufficient strength to deviated solar 
wind particles into the detector. In addition, there is still 
no known acceleration mechanism that could produce the 3500 eV 
particles that are observed. There is no reason to suspect 
that the spectrum ends abruptly at 3500 eV as well. Counter¬ 
streaming protons of energy 3 to 6 keV have been observed in 
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this region of space by other investigators, viz. Asbridge 
et al. [1968]. Lin et al., [1974] have observed counter- 
streaming protons in the energy range 30 to 100 keV. Although 
solar wind protons in the energy range 5 < £ <50 keV known 
as suprathermal solar wind protons have been observed [Frank, 
1970], these energetic proton events are not characteristic of 
the quiet solar wind. Since these suprathermal protons come 
from the solar direction one would again have to invoke some 
mechanism to deviate the particles into the detector. 

c.) cosmic rays-Very high energy protons ( £ > 50 MeV) 

can penetrate the case of the SIDE instrument and activate 
the channeltron causing counts to appear in whatever energy 
channel is being sampled at the time. The cosmic ray flux 
in the interplanetary medium is responsible for part of the 
background noise during instrument daytime and is responsible 
for most of the background during instrument nightime. These 
counts in random channels are generally very low (1 or 2), and 
are only enhanced during unusual solar proton events such as 
solar flares as discussed by Medrano [1973]. That is, cosmic 
ray protons produce a random background of very low intensity 
in quiet times and ir disturbed times the background counts 
rise independent of energy channel making their identification 


easy. 
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d.) other sources-At the time of this writing no lunar 
bow shock capable of producing 3 keV protons has been observ¬ 
ed. Other complicated solar wind-lunar interactions may occur 
to produce these particles, but such processes have not been 
reported. 

The above considerations alone indicate strongly that the 
high energy particles observed in this region are indeed 
protons from the bow shock since it is the most likely other 
source. Other corroborating evidence that the high energy 
particles are from the bow shock can be summarized as follows. 

a. ) Protons of the same energies as the observed ones have 
been previously observed in this region by several investiga¬ 
tors so it is not unreasonable to expect that the SIDE would 
detect them. The instrument looks back at the shock in this 
region of the orbit. 

b. ) Results of the data analysis of this thesis indicate 
that these high energy particles are influenced by the inter¬ 
planetary magnetic field in such a way as to be consistent 
with the proposed EXB drift model. 

It is apparent that bow shock particles are a dominant 
feature of the data observed in the bow shock regions which 
suggests a reclassification of the recurring ion events 
discussed by Lindeman [1971,1973]. Lindeman classified the 
ion spectra in the extra-magnetospheric regions of the lunar 
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orbit into five basic types, two of which are clearly ident¬ 
ifiable with bow shock protons. The Type I spectrum has a 
single peak from 2000 eV to 3500 eV and Type II has a double 
peak, the first peak from 30 eV to 250 eV and the second 
from 2000 eV to 3500 eV, (see figure 6-4). S ub sequent inves¬ 
tigation by Lindeman [1973], has shown that the low energy e- 
vents in the bow shock region can be classified into three 
types. 

a. ) VXB events 10 eV to 500 eV 

b. ) Lunar surface potenial (LSP) events 10 eV to 100 eV 

c. ) Instrument resonance events (IRE's) £.50 eV 

In light of these findings, the low energy peak of the 

Type II event is easily identified as being due to LSP and 
VXB events near the terminator with the high energy peak be¬ 
ing explained by the presence of bow shock protons in this 
region. The Type II event is merely a Type I event with the 
contribution to the spectrum of the LSP and VXB events. 

The ion data observed in the bow shock region can there¬ 
fore be classified as consisting of four basic types of events 
a.,b.and c. above and 

d. ) Bow shock proton (BSP) events £ a 750 e V 

Based on the present research it is possible then to state 
a general criterion for identifying BSP events in the data. 
During passage of the moon through the bow shock region, any 
observed particles having energy greater than or about equal 



Ficrure 6-4. Recurring ion events Types I and II. Type I 
has only a high energy peak while Type II has a low energy 
peak in addition to the high energy peak, [redrawn after 
Lindeman (1971)] 
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to 750 eV are probably bow shock protons, with the exception 
of clearly identifiable high energy events associated with 
unusual activity in the interplanetary medium. 

Figure 6-5 shows representative counting rate spectra in 
the bow shock region near the sunset terminator with two of 
the above types of events easily identifiable. LSP's and BSP r s 
are readily apparent and a tentative identification of V X B 
events can be made. The marginality in the identification 
of the V X B events could be easily removed by looking at the 
interplanetary field behavior for this period and remembering 
that B out of the ecliptic implies V X B events and B in the 
ecliptic implies BSP events. 
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Figure 6-5. Representative part of the bow shock region 
showing varied ion spectra near the sunset terminator. 
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6.2 Further Research 

Limitations on the accuracy of the present analysis are 
due primarily to a lack of sufficient simultaneous magnetic 
field and ion data. An analysis of the bow shock regions 
carried out over several lunations would much enhance the 
significance of the results and would likely improve the agree¬ 
ment of the data with what is expected from the model. This 
can be seen by noting the better correlation achieved in the 
outbound region/ (see figure 5-5) f where three days of data 
were analyzed as opposed to two days on the inbound side. 

The amount of simultaneous particle and field data available 
in any one lunation is also limited by the data gaps in the 
field data incurred while Explorer 35 is behind the moon. 

The lack of an onboard recorder causes data gaps of about 5 
hours on each orbit of the satellite. 

Other problems exist such as the unknown spatial distri¬ 
bution of the field, an important factor in determining the 
resultant guiding center trajectory of a particle. Not much 
can be done about this problem but it appears that the var¬ 
iation of the field over a distance comparable to the proton 
gyroradius is not significant since particles do seem to make 
it to the moon via the EXB mechanism. Likewise the problem 
of including the effects of the local surface fields is very 
difficult in that the spatial configuration and extent is not 
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well known, but again these fields dp not seem to play a 
significant role from the self consistent point of view. That 
is, if one assumes a model and observes the behavior predict¬ 
ed by the model then other influences must not be too signifi¬ 
cant. 

In some instances the field is also highly time varying 
and no account of this has been taken in the present analysis. 
If the field is varying rapidly with respect to a character¬ 
istic time of the motion, say the gyroperiod, ( 12 sec) , the 

guiding center treatment breaks down. There were several in¬ 
stances in the data where the field was stable with respect 
to a gyroperiod but the inclusion of times where the field 
was varying rapidly would detract from the accuracy of the 
results. An obvious way to improve the results would be to 
confine the analysis to times of stable field conditions. 

As noted above, the bow shock reyr'ons are of a highly 
diverse character in ion spectra with bow shock protons con¬ 
tributing a large portion of the observed ions but certainly 
not all of them. This would imply that any correlation of 
particles observed with field configuration would be further 
obscured by including particles whose observability did not 
depend on the field configuration. All events were included 
in this analysis. Obviously, the contribution by the LSP's, 
the IRE's and any other event associated with the lunar atmos- 
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phere would be detrimental to the results. This point can 
be clearly seen by noting the effect of including all of 
the ion events on the distributions shown in figures 5-5 and 
5-3. In particular, the inclusion of the VXB events inflates 
the upper end of the distribution in the outbound region and 
the lower end of the distribution in the inbound region. For 
this particular type of event the contributions to the dist¬ 
ributions cannot really be considered detrimental since VXB 
events do depend on the field configuration and are easily 
distinguishable from BSP events. 

Again by considering isolated events where the field is 
stable and where particles from other known sources can be 
removed from the analysis, e.g. near the terminator, the cor¬ 
relation would be greatly improved. 

With the availability of good data from all three SIDE'S 
and Explorer 35, the directional nature of the proposed model 
could be more effectively tested. The wide variation in long¬ 
itude of the look direction of the three detectors would al¬ 
low a cross correlation in dependence on longitude of the ob¬ 
servability of ions. The use of three instruments in this 
manner could also be helpful in determining the effective 
size of the look cone as discussed in chapter 5. 

Any investigation such as the one done for this thesis 

yields many new and interesting areas for further 
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work. Two of the related areas of special interest are, the 
possible acceleration mechanisms associated with the bow shock 
itself and a study of the nature of the source points on the 
shock surface. If a very good fit to the real data can be 
achieved using a particular type of source function, then one 
is led to suspect that that source function probably approxi¬ 
mates the real one. Knowing the source function at the shock 
surface could be useful in determining which if any of various 
models for the acceleration mechanism is correct. Two of the 
major models for the acceleration mechanism are reviewed in 
the Appendices. Specular reflection at the shock and subse¬ 
quent energization by the interplanetary electric field has 
been proposed by Sonnerup [1969] and the mechanism of Fermi 
acceleration has been developed by Jokipii [1966]. 

However, the primary mission of the SIDE is to gather 
data that can be used to learn about the lunar atmosphere 
which suggests that the real value of this research and any 
further research in this area is that it aids in understand¬ 
ing the complicated spectrum of events occurring in the bow 
shock regions of the lunar orbit. These bow shock regions 
are the most fertile in the data for providing information 
on the lunar atmosphere and certainly warrant further inves¬ 
tigation . 
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Chapter 7 


7.1 Summary 

The basic objectives of the research done for this thesis 
were outlined in chapter 1 and will be briefly restated here. 

a. ) to determine if the particles of energy >750 eV ob¬ 
served in the pre and post bow shock regions are indeed bow 
shock particles. 

b. ) to establish a criterion for the identification of 
these particles. 

c. ) to propose and develop a physical model to explain 
the presence of the bow shock protons if they are found in 
the data. 

Such a model based on EXB drift trajectories was developed 
in chapter 3 leading to the following basic conclusions about 
the behavior of bow shock protons in this region. 

a. ) the observability of bow shock particles at the moon 
by this mechanism is highly dependent on the configuration 
of the interplanetary magnetic field. 

b. ) there is an inherent dusk-dawn asymmetry in the 
model that should be observable in the data. 

It was also pointed out that the present model/ although 
similar to that proposed by Asbridge at al. [1968], is a more 
realistic one in that it describes all particles leaving the 
shock/ not just those with a special velocity and pitch angle. 
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The details of the computer model and its connection to 
the physical model were discussed in chapter 4 where criteria 
for the observability of bow shock particles at the moon were 
established. For a particle to be observed at the moon by the 
EXB mechanism, it must, 

a. ) come from the right direction and 

b. ) have a source point on the shock surface. 

The plausibility of the physical model of chapter 3 was 
established by the data analysis, the results of which are 
displayed in chapter 5. The major conclusions drawn from 
interpretation of the data are, 

a. ) analysis, independent of the model, shows that the 
data behave as one would expect from the model. 

b. ) synthetic data generated by the computer model agree 
semi-quantitatively with the observed data. 

c. ) behavior of the computer model coupled with consider¬ 
ations of actual particle trajectories indicates that the 
effective size of the look cone of the detector is substanti¬ 
ally larger than the nominal value when applied to bow shock 
protons. 

Arguments against other possible sources for the observ¬ 
ed ions of energy >750 eV are presented in chapter 6. 

These other sources are, 

a.) lunar atmospheric ions 
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b. ) solar wind ions 

c. ) cosmic rays 

d. ) lunar bow shock or other unknown sources 

and were rejected primarily on energy considerations. Lunar 
atmospheric ions and solar wind ions have much lower ener¬ 
gies in general than the observed moderate to high energy 
ions and it was pointed out that no known acceleration mech¬ 
anism at the moon is capable of producing particles of the 
observed energies. Cosmic rays were dismissed based on the 
well established character of penetrating particle events 
which are unlike bow shock proton spectra. 

The value of isolating and identifying bow shock protons 
in the data was also discussed in chapter 6 where it was noted 
that the usefulness of the bow shock region for studying the 
lunar atmosphere can be enhanced by separating BSP events from 
the lunar atmosphere events. A general criterion for ident¬ 
ifying bow shock protons was given; particles of energy great¬ 
er than or about equal to 750 eV and found in the pre and 
post bow shock crossing are probably associated with the 
bow shock. 

7.2 Conclusions 

a. ) The protons observed in the bow shock regions with 
energies > 7SO ev are from the bow shock and can readily 
be identified in the data. 

b. ) The proposed EXB drift model is a viable one as long 
as its limitations are kept in mind. 
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Appendix 1 

Al.1 Define an orthogonal coordinate system S and an orth¬ 
ogonal coordinate system S’ that has the same origin as S 
but is oriented arbitrarily with respect to S otherwise. It 
is possible to make the axes of system S coincident with 
system S' by three sucessive counterclockwise rotations of 
the S axes. The angles through which one must rotate to per¬ 
form such a "transformation" of the coordinate axes are known 
as Euler angles. The angles <£ , i}j , and 0 of figure Al.l are 
the Euler angles for the transformation from system S to S 
and the transformation matrix and its inverse are shown on 
the next page. 

Al.2 Notes concerning chapter 4. 

In transforming from system S’' to S sq the matrix of 
transformation was ( SBA ) and the position of the moon in the 
SE system was found to be 

"r = SBA ~r 1 ' Al. 3 

se 

The vector ~r 11 has only an x-component in system S 1 ' 
since the moon is on the x''-axis at a distance r from the 
origin. That is, 


r 1 ' = r x ' ' 

or represented as a column matn 


r" = 0 


Al. 4 


0 




138 


8 

•3 8 


u u -e- 

<33 <35 m 

01 01 O 

o o u 

° ° <35 

+ + R 

*©- *0- tn 

r a i 

■rl 


ffl III 

O O R 

CJ o -H 



139 



Figure Al.l Euler angles 6 , ^ and <P . The coordi¬ 

nate transformations between the primed and the unpriraed 
systems can be conveniently expressed in terms of the 
Euler angles. (Goldstein, 1950) 





2 



Figure Al.l 


Figure Al.l 
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This column matrix, when multiplied by ( SBA) from the left 
gives 




A1.7 


Hence equations 4-8 
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Appendix 2 

Derivation of equation s 3-2 through 3—9 . . 

Consider a charged particle moving in a region in which 
there exists electric and magnetic fields. The force on the 
particle due to the electric field is, 

F e = q E A2.1 

and the force due to the magnetic field is 

~F b = q V X ~B A2.2 

so that the equation of motion (eq. 3-1) is 


m-ir" = q ( "E + *V X "b ) 
at 


or in component form* 


dV 
—-x 

= q 

( E 

+ 

V B 

- V B ) 

A2.4 

dt 


X 


y z 

z y 


dv 

d£ y 

= <1 

1 E y 

+ 

V B 

z X 

- V B ) 
x z 

A2.5 

dV 
- 7 ~z 

= q 

( E 

+ 

V B 

- V B ) 

A2.6 

dt 


z 


x y 

y x y 



Since the electric field is given by 


E is always perpendicular to B. So without loss of generality 




we can take B=Ez and "E=Ey. Equations 4,5, and 6 of this- 
section then become. 


=a - V y B 


A2.8 


= q E - q V B A2/9 


= 0 


Taking the time derivative of equation A2.8 and substituting 
it into equation A2.9 gives 


m(m/qB) -^x = q E - q V B A2.ll 

at 2 x 


2 

——^x + (qB/m) 2 V =(q 2 E B/ra 2 ) A2.12 


Defining 


fi) = q B / m 


and 


2 2 
K = (q E B/ m ) 


equation A2.12 becomes 
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Equation A2.14 is a well known second order, first degree 
linear differential equation with constant coefficients whose 
solution is of the form 

V x (t) = A cos wt + D sin cot + C A2.15 

The constant C can be determined by substituting equation A2.15 
into equation A2.14. i.e. 

V = - A sin cot + 

x 2 

V = - co Acos tut 

y 2 

- - (V x - C) 

so that 

\ + *> 2 \ - K A2.17 

becomes 

C u 2 = K A2.18 

or 

C = K/o) 2 = E / B A2.19 

A and D can be determined by applying the boundary conditions 
at t=0 and equation A2.8. 

A = V xq - E / B A2.20 

and from equation A2.16 and A2.8 


D cos a t 
2 ' 

co Dsincot A2.16 
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'V (0) = <B D 3 qB/m V 
X yO 


A2.21 


or 


A2.22 


and hence. 


V x (t) - (V xo - E/B) coscot + V yo sin cat + E/B A2.23 
Similarly, 

V y (t) = V yQ cos cot - (V xq - E/B) sin cot A2.24 


and from eq. $.2.10, 


V (t) = V 


A2.25 


The quantity E/B wiil be shown to be the magnitude of 
the drift velocity. Therefore the equations for the velocity 
are 

V (t) = (V - V) coscot + V sin cot + v. 

x xo d yo d 

V (t) = v cos cot - (V - V,) sin cot 

y yo xo d 

V (t) - V 

Z 20 

Integrating eqs. 26 through 28 of this section and ap¬ 
plying the condition that the particle start cut at the ori¬ 
gin yields, 


A2.26 

A2.27 

A2.28 




145 


X(t) 

Y(t) 

Z (t) 

Equations 29 and 30 are the parametric equations of a cycloid 
in the x-y plane arid the z-motion is independent of the x-y 
motion. The transformation from a fixed coordinate system to 
one that is moving parallel to the x-axit, with velocity V p is 

X' = X - V^t 
d 

Y’ = Y A2.32 

Z* = Z 


- j^(V xo -V d ) /cu^ sin < 0 1 + (V^/co ) (1-cos cu t) + 


t A2.29 


= (V^/cu) sin w t + 


Q -V d ) /Jj (cos to 


t - 1) A2.30 

A2.31 


This implies that eqs. A2.30 and A2.31 are unchanged in 
such a transformation but equation A2.29 becomes/ 


X'(t) = 


- ( (v x 0 - v a>^)‘ 


in cut + V yo /&> (1-cos cut) 


Let A = (V -V )/cu and D = V /co 
xo d yo 


p=<ut then/ 


X* = A sin<£ + D (1 - cos<£) 
Y' = D sin<£ + A (cos <p -1) 


Rewriting eqs. A2.33, 
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(XV - D) - A sin<£ - D cos <p 
(Y 1 -(-A)) . — D sin<£ + A cos 9 ! 


Defining the angle, 

Q = tan" 1 (A/D) A2.35 

then equations A2.34 become, 


(X’ - D) = ^Ja 2 + d 2 cos ( Q + <f> ) 
l ^fi 2 + D 2 , 


(Y' -(-A)) ~=\JA 2 + D 2 sin( 9 + <f> ) 


(X ' - D ) 2 + (Y 1 - (-A) ) 2 = R 2 


Apparently the equations A2.29 and A2.30 represent cir¬ 
cular motion when transformed to a system that is moving par¬ 
allel to the x-axis at velocity V^. The center of the circle 
in the moving frame is 


X = D = V /o» 
C y O 
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and it has radius, 


R = 1/(0 ( (V -V,) 2 + V 2 ) 1/2 = mVfV qB A2.39 

xo a yo 


where is the perpendicular velocity in the moving frame. 
Hence equations 3-5 and 3-6. 

The particle appears to move under the influence of a 
magnetic field only in the moving frame. That is, we have 
transformed into the frame in which there is no electric field. 

= y( |+VX B ) A2.40 

Consider the relation, 

B X (4.X B) = B 2 ^- (B'E|} B A2.41 

E ± = - (E x X 13}XB/B 2 
but from eq. A2.40, 


A2 0 42 


Evidently, 


E. = - V X B 


V = (E X “1) / B 


v/hich is equivalent to equation 3-8. 
The magnitude of V is 
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Hence the particle obeying eg. A2.3 seeks circular motion 
in a frame in which the electric field is zero, i.e. one 
moving with velocity v such that the component of the vel¬ 
ocity that is perpendicular to the field is equal to the 
drift velocity. Obviously, the particle can move in the z- 
direction as well as in the x-y plane, but any velocity it 
has in the z-direction is velocity parallel to the field and 
it will still be moving in a frame in which the electric 
field is zero. Again, since the z-motion and the x-y planar 
motion are uncoupled, the total motion can be described as 
cycloidal in the x-y plane and uniform translation of that 
motion in the z-direction. The resultant velocity vector of 
the particle's' guiding center is 

*V = V z + V, x A2.46 

r zo d 

which is equivalent to eq. 3-7. Clearly, the particle moving 
in a frame with velocity^, sees no electric field. 

- 7 < \+ (V 2 Q Z+ V d x) X B z) 

= y ( Ej_+ V d B x X z) A2.47 

= r ( %.-%) - 0 

Since the perpendicular E field is due to the solar wind 

velocity "v in the present treatment, using eq. A2.44 and 
sw 

'E=-\ij 3w X B it is easy to see that 
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"v = V - (V * B) B A2.48 

d sw sw 

_^ ^ ^ -A 2 

v;= ( -(V X B) X B ) / B A2.49 

d sw 

or using the identity , 

A X(B X C)= (B*A) C - (C-A) B A2.50 

% = <^sw E ' - t 5 / 

% = ***- A2 ’ 51 


hence equation 3-9. 
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Appendix 3 

& 3 . i a representative correlation program is discussed in 
this section, viz. the program that produced figure 5-6. The 
basic idea behind all of the correlation programs was discus¬ 
sed in chapter 5 and will not be repeated here. Only the de¬ 
tails of the procedure used to obtain the results are relevant 
to this discussion. 

Program A3.1 does the following: 

a. ) reads the Explorer 35 magnetometer tape to find the 
start time 

b. ) reads the SIDE data tape to find start time 

c. ) reads a record of Explorer 35 data to determine the 
value of the field latitude for that 82 second interval and 
decides the appropriate bin number for that latitude (i.e. for 
the 180 values of field latitude from -90° to +90°, 90 bins 

of 2° each are used and are numbered from 1 to 90). 

d. ) after selecting the bin number, the average integral 
flux over the same 82 second interval is calculated from the 
SIDE data tape and this number is deposited in memory at that 
bin location. 

e. ) The processes c. and d. above are repeated until the 
period of interest has been evaluated. Each time flux is ad¬ 
ded to a particular bin, a counter for that bin is incremented 
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so that the number of times that bin was encountered is 
known. 

f. ) at the'end of the period of interest, the total 
accumulated flux in each bin is divided by the number of 
times that that bin was encountered. Thus the number ob¬ 
tained in that bin represents the average integral flux 
observed at a given latitude over the period of interest. 

g. ) The results of all 90 bins are then dumped onto 
magnetic tape and are subsequently plotted to obtain figure 
5-6. 


The details of the calculation of the integral flux are 
as follows. Figure A3.1 shows the meaning of the parameters 
relevant to the following discussion and the approximate na¬ 
ture of the integration over the energy spectrum. The appro¬ 
priate values of the parameters for the Apollo 14 SIDE are 
given in table A3.1. 

For the ith energy channel, the number of counts C* will 

be given by . 

e = R At R , = count rate (sec ) 

i i i . 

At = 1.13 sec 

the count rate is related to the differential energy flux 

by ■ ■ . . . 
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Figure A3JL Illustration of the technique used to cal¬ 
culate integral flux. The number that represents the 
integral flux is given by summing the areas of the unshad¬ 
ed rectangles defined by AE^ and (dj/dE). 




Figure A3. 
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50.0 
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3.2 

76.0 
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9.0 

140.0 
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4.7 

230.0 

237.5 

7.9 

180.0 

225.0 

9.8 

190,0 

225.0 

3.3 

250.0 

262.5 

9.9 

320.0 

287.5 

8.8 
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262.5 

7.1 
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Table A3. 1 

peak response 
7 


1 

10 

2 

20 

3 

3 b 

4 

50 

5 

70 

6 

100 

7 

250 

8 

500 

9 

750 

10 

1000 

11 

1250 

12 

1500 

13 

1750 

14 

2000 

15 

2250 

16 

2500 

17 

2750 

18 

3000 

19 

3250 

20 

3500 
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where/ 2 

G i = G i^ (cm " str ^ 

Ae^ = energy passband (FWHM) 

The flux, integrated over the energy spectrum is j. 



AE Wi = width of the i th integration interval 


or in terms of the known parameters, 

20 / , ■> 
J =y'fl^5 ! i - 1/2 C 1 AE “1 + °2Q AE «20 

" At AE i G i 1 At AE 1 G 1 At AE 20 G 20> 

(See figure A3.1 to understand the second term) 


The values of Ae^. were determined by taking half of the 


energy difference between the i th peak and the i-tl th peak 



1/2 K+i 
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A3.2 Program A3.2 is the program used to generate the syn¬ 
thetic data of chapter 5. The theory underlying the gener¬ 
ation of these data was covered in chapter 4. The program 
presented here is the SDS 910 version. The IBM version is 
somewhat lengthy and will not be included. It is basically 
the same as the 910 version except that the ephemeris data 
are calculated in that program. However, the tests relevant 
to the acceptance criteria are the same. 

To facilitate the understanding of program A3.2, it should 
suffice to define the parameters relevant to the theory of 
chapter 4. 

T=kinetic energy of the particle 

D=d=distance from the center of the earth to the nose 
of the shock. 

CAPD=D=distance from the center of the earth to the wing 
of the shock measured along the y-axis. 

CR=r= s position vector of the moon 

APX,APY,APZ, components of the detector look vector in 
a solar ecliptic system d 

VR=V r =resultant velocity vector for the guiding center 

VDOTAP= V • D 
r 

B2 and AC4 represent equations 4-22 


HGT=counts 


A3.1 
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DIMENSION IRY! 2&) * BIN ( 91 ),[T.IFE(5) ,IDIV!91>,FLX(10). IW120) 

CALL lTY PE(lTIME ) 

C TYPE IN START TIME 

CALL ITYPEURL'N) 

C TYPE IN RUNTIME IN HOURS 
CALL lTYPE (I M) 

C TYPE IN MODIFIED GEOMETRIC FACTORS 
HTODAY=1. /24. 

TODAY^HTODAY/60. 

ST0DAY=HT0DAY/36 
DGAP=i . *TODAY 

T8=FLtJAT( niMEa 1 ) > + F.LC ATU I TIME ! 3 > ) *FTCDAY+FLOAT! [TIME!4) >*TT3DAY+FL 
1 OA T( IT [ME! 5) )*ST OOAY 
TEND=TB+RUNTLK 
DO 1 K ” 1,9 1 
dl N! K J =0 . 

1 ICIV!K)=C 

ASSIGN 2 TO [PAR 
ASSIGN 3 TO IE OF 
l CONTINUE 

CALL SETUP! 1. 28, IRY. IHECI, IPAR.IEOF ) 

C FIND START TUE ON EXPLORER JS TiAFE 


>2 ; 


READ INPUT TAPE 2,16, ID A. Y,[HR, II" l N , ISEC 
FORMAT (4 [3) 

IDAY-1DAY+1 

IF! I DAY-IT [ME! 2) II 7, 18.2 i 
IF(I rfi-ITIMElJ)) 17,19.21 
IF!IMIN-ITIME<4> )17,2C.2 l 
IF (ISEC—IT I ME! 5) 117, 21,2 l 

T8=FL0AT( IDAY)+FLOAT< I hr MHT.OLAY-FLGAT! IM[N)*T OCAY+FL OAT ! I SEC ) *5 TL 
1DAY 

FIND START TIME ON SIDE TAPE 
CALL I NPUT 

T= FLOAT I IRY! 2 H +FLOAT l IR Y!3))*HTOD A Y+FL OAT ( [RY14U*T QD AY AFLOAT ! [R 


lYi5 > »/ l n JO. ) fSTODAY 
25 IF !T-T81 23.2 4, 24 

2 3 T-T+24.16*ST ODAY 

CALL INPUT 
IF!IRY(Z>)25,22,22 
24 TPLUS=T + 24.16*ST GD AY 

TL AST= TB 

C 8CTH TAPES ARE ON START TIME HERE 
C READ EXP 35 TAPE 

1- READ INPUT TAPE 2,2 6 , 10 A Y . I HR. IM IN . 


OF POOR QUALM 
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26 IFCRHAKAIJ.32X.F6-1} __ 

TB=FLOATII-FLOAT I1 HR I •HOQtAY.+ F LOATt IHlNI *T 00AT-FLOAT t ISEC) ''STC 
1DAY 

gap=tb-tlast 
IFtGAP >30, 30.59 
5C IF(GAP-OGAP) 48.4 8, 51 
C CALCULATE BIN NUMBER 
48 I=BL AT/2. 

1=1+46 

N=Q 

SUH=0. 

C CALCULATE INTEGRAL FLUX 
29 CALL INPUT 

N=N+1 
CT-O. 

DO 4 J = 1.20 
IFtlRYtJ+5>)6 .4.4 

/* rr-rT+FI OATt IRY ( J+ 5) ) * FL DAT tlW(21“JlL/10- _ 

CT =CT-( FLOATt IRY t6> > *1 Ml 2 . >+ FLOAT( IP Yt 25 ) > *IM< 1 > >/ 20- 
62 SliM-SUM + CT 

IF(IRY < 2) >27.28,28 
27 TP=TPL US 

„ yp=FLDAT ( JR YI 21 > *FL 0 AT ( IRTI 31 I*HTO DAY-FLOAT! IRY! 41 I —TOD AY-^ LOAT! I 

IRY(5)>/1- - 0.)*ST UDAY 

}3 TPLUS=TP+24. L6*S TODAY 

IF(TP-TB)29,31, 3 L 

31 SUH=SUH/FLOAT (N) 

BINt n=BINtl )+SUM 
ID IV < l >=ID IV l l > + 1 
TLAST=TB 

[F(TB-TEND)30,32.32 

32 CONTINUE 

CALL I TYPE(I MATE > 

C MCUNT BLANK TAPE ON UNIT 1 
ID C 13 K= 1,81.13 
DO 5 7 J = l. 1 • 

I =K+ J- 1 

IF (TO I V ( I >) 14.14., 56 
14 FLXt J> =C . 

DO TO 57 

56 FLXCJ>=BIN(I >/FLOAT( ICIVKI)} 

57 ICCNTINUE ^ 

13 . WRITE OUTPUT TAPE 1.58,1 ICIV tw ). J-K.K xj 

5 8 FORMAT (8X, 10 11 A > 

END FILE 1 
3 ICCNTINUE 

END 





A3.2 


158 


C PROGRAM A3..? 

C J. BENSON RICE UNIVERSITY- 
C IP (1 > = D I ST ANCE TO NOSE OF SHOCK iIA RE 

C IP (2 1/IP< 3 ) = DISTANCE TO MOON- AT BOW SHOCK CROSSING IN RE 
C IPU>/HPt&.>*SCLAR ECL LONG AT BOH SHOCK CROSSING IN RADIANS 
C II F( 6 ) = SOLAR WIND VELOCITY INI KM/SEC 
C IP(7 ) = SIZE OF THE LOOK CONE HU DEGREES 
C IP(8) = SHAD0W ANGLE IN DEGREES 

C ’.SHADOW ANGLE IS USED. TO IAV01D UNECESSARY TESTING 
C IPC9 1= NUMBER OF RECORDS iTC BE PROCESSED AT A TIME 
C IF C1 I =SOURC E NUMBER 

DIMENSION T( l .IP t 10), IRi24 ) »R (12 ), LCUTt 28) 

EQUIVALENCE(IR C 2J,R< UI 
ASSIGN 40 TO IPAR 
ASSIGN 4 1 TO I EOF 

C CALCULATE 2“ VALUES OF T-HE ENERGY 
IDC 2 K - I t L 4 

2 TC21-K)=350C.-CK-i)*25G. 

T( 6) = 1 * 

IDC 3 K = l,3 

3 T ( 6—K) =7 .-( K— 1J *20. 

U 11=1'". 

T C 2 > = 2 •' . 

CALU ITYPE(IP) 

C TYPE IN INPUT PARAMETERS 
O-FLOATC IP( 1 II 

RE = FLnAT(IP(2)l/FL0AT(IPC :IJ 
SELAHS =FLQAT(IP(4 > >/FLOATlIPC 5)) 

VS=FL0AT IIFC &) ) 

SLA=FLOAT(IPC 7 )) 

SHAD=FLOAT(IP<8)) 

BETA=. 8726 
PI =3.14159 
PI2=P1/2. 

'TCRAD-P 1/180. 

SL A= SL A*TCHAD 
CLA=SlNCPI 2-SL A) 

S6S s SINt SELA8S) 

CBS = SIN< PI2-SELA8S) 

C CALCULATE VALUE OF WING DISTANCE CAPD 

CAPD=SOR T( C-0*RE*REv SBS*.SBS) /<RE*C BS-D)) 

SB-S IN CBETA) 

C8=S IN CP I2-B ETA) 

4" CALL SETUP! 1 • 24, III, IR6C. IPAR .1 EOF) 

IDC 7 1 = i , fP( 

C READ EPHEMERI5 DATA 
CALL INPUT 
7 CONTINUE 
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C TEST SHADOW ANGLE 

IFtR (1 D-SHADJ1* l. 5 
1 'DO 6 J = 1»20 
6 iGUT(J +5 ) = 0 

GC TO 4 

5 IF(R 1111-180.>9, ib,16 

16 IF(R(1 I)-SHAD-I8u.>1.1.9 
9 ICCNTINUE 

C CALCULATE SOLAR ECLIPTIC ANGULAR! PARAMETERS 
SE LD=S IN<R(7 ) ) 

CELD =5 IN(PI2-RC7)> 
sseloc=sin(R tan 
CS EL QD = S IN (P 12-R t 8 )) 

SE LA =S IN t R ( 9 )) 

CEL A=SIN(PI2-R(9 1 ) 

SSELOA=SIMR (12) I 
CSEL0A = SIN(PI2-R(l r ') > 

SBLQN=SIN(R( 11)* TORA D) 

C6L0N= SI N( PI 2-R( li)*TCRAD) 

SB LA r ='S INC R ( .2)* TOR A 0) 

CBLAT = SIN(PI 2-R( i2)*TCRAD) 

BX =C8L AT *C BL ON 
BY=CBLAT*SBLUN 
BZ = S BL AT 
IF (8 Y ) 17,19, 19 
19 BX=-BX 
8Y = -BY 
BZ =~ BZ 

C CALCULATE R 

17 CR = SQRTCR(4 ) *R( 4 )+R( S ) *R t :) +R *R < 11 ) 

CR X = CR *C EL A* CS EL UA 
CKY = CR *C EL A* SSEL UA 
CR Z=CH*SELA 

C CALCULATE DETECTOR L:UCK \l ECT'Of. 

APX=-C EL D*CS EL OD 
APY=-CELD*SSEL OD 
APZ=-SELD 

DlfNC = VS*CE*C6LAT *cbugi\ - v s *sb *CBL at *SBLON 
C EVALUATE 2? ENERGIES 
DO 4 K=l,20 

IVF = SOKT< T(K ) * 1.6 2 E- 22 ) * i."9 it £1 2 
HGT = .. 

C EVALUATE if PITCH ANGLES 
DO 13 J-1,10 
PA = < J- *iO. * TOR AO 
SP A = SIN(PA) 

CPA=SQKT< l SPA* SPA) 

VP =VP* CP A 


PM6MAU PAGE IS 
lOl.POOB (QUALITY 
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C CALCULATE RESULTANT VELOCITY! VEC.TCR 
VRX = -VS*CB + BX*DUNG*BX*VP 
VRY-VS *SB+BY *DUNG+8Y *VP 
VRZ=BZ*(DLNG+VP) 

VR=SQRT( VRX* VRX+ VRY* VR Y+ VRZ* VRZ>. 

C CALCULATE ANGLE BETWEEN VR ANC DEI LOO 
VDOT AP=( VR X* APX+VRY * AP Y+VRZ * APZ ) / V 
IF (VliOTAP-CL A) 13.10. K 
io iFtVDOTAP-l. )18. i3.1 3 

C TEST FOR SOURCE POINT 

18 CAPA=.+<<VRZ*VRZ>/tVRY*VRY>) 

CAPB=( <C APD*CAPD »/D)*< VRX/VRY1-2* < 


K VEC AND TEST 
R 


CP'Y*» VRZ/VRY >*< VRZ/VRY )-CRZ♦ C VRZ 


1 < PARF =*t (C R Y*V RZ)/ VR Y ) *( (CRY*VRZ) /VR Y) 

BELCH=-CRY * ( < Z*CRZ*VRZ) /IWRY+C CAPO* CAPO* VBX »/ID* VRY1) 

ERUCT= <C APO* CAPD *CRX ) / D+ CRZ* CRZ-CA PD *CAPD 
CA PC =B AR F+BELCH+ERUCT 
B2 = C AP B* CAPB 
AC4=^*CAPA*CAPC 
IF (8 2- AC 4) 13 . it. il 
i t [F (PA) 13. I 4, 15 
C FILL CHANNELS WITH APPROPRIATE COLNTS 
HGT=‘ IP ( U>*1 .*T OR AD 
GO TO 13 

HGT“HGT+IP< U I *S INIPAI *1 , .*TORAC 
ICCNTINUE 
10 UT ( 2 6- K) =H GT 
CONT INUE 


11 


15 


IOUTCi)=-l 
11 C UT t 2 > = R t 1) 

IOUTn ) = R( 2) 

IQUT(4 >=R(3) 

11 CUT(5)= j 
iDC 8 J = 1 • 3 

8 CALl'sETUPC2.28.ICUT,I REC.IPAR.IEOF ) 

C WRITE A SPECTRUM ON MAG TAPE 
CALL OUTPUT 
GO TU <*'■ 

41 CONTINUE 

END 
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Appendix 4 


Brief review of the acceleration mechanisms. 

A4.1 Sonnerup model - The basic idea of the Sonnerup 
model is simple specular reflection from the shock and sub¬ 
sequent energization of the particle by the interplanetary 
electric field. To get an order of magnitude estimate of the 
increase in energy after reflection we can calculate the ra¬ 
tio of the reflected energy to the initial energy. Assuming 
the reflected energy is given by the initial energy plus the 
energy given to the particle by the electric field over one 
half of a cyclotron orbit. 


f = 1/2 mV 2 + q E nr 
r sw M IP 

where 

r - mV sw /qB 

and 

= V B sin 
IP sw 

^ = angle between V and B 
sw 

2 2 

e = 1/2 mV + 77 mV si nifr 

r sw sw 

so that. 


A4.1 


A4.2 


A4.4 


2 2 
(mV 77 sin iA ) / (1/2 mV ) 


A4.5 
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e r / - 1 + 27rsimA 
for a typical value of (\fj = v/ 4 ) t 


A4.6 


A4.7 


That is, a particle could gain roughly 5.5 times its incoming 
energy by this process. For a typical bulk flow velocity of 
about 400 km/sec, the reflected particle energy would be, by 
eq. A4.7, about 4.7 keV. This is the right order of magnitude 
for the observed particles in the upstream region. For dif¬ 
ferent values of the energy gained varies from 0 to about 
7.28 times the initial energy. 

Sonnerup treats this problem in a more rigorous fashion 
in his 1969 paper where he considers that the particle energy 
is conserved in the frame of reference in which the electric 
field vanishes. The results of this more rigorous treatment 
will be stated here. By considering the reflection process in 
the frame of zero E field and the rest frame of the shock, 
Sonnerup obtains, 

/ _ ^ (1+5) cos ( <£ + y ) sin(?A - 4>) sin( ^ + y) _ . „ 

V e i “ 1 + 2 -(l+fi)cos 2 (*-<£) A4 * 8 

where the angles are those shown in figure A4.1 and the def¬ 
initions of /x and S are as follows. 


V 1 i /v i 


s= v »i 


A4.9 
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Figure A4.1a . Vector diagram showing how the incident 
particle guiding center velocity vector, , is decom¬ 
posed into a component, *V J( i , along the field and a com¬ 
ponent V t along the shock front. 

Figure A4.1b. Similarly the velocity vector of the re¬ 
flected particle is decomposed into V J( r along the field 
and V t along the shock front. The angle 'P between the 
field direction and the solar wind vector V gw is the 
garden hose angle. The angle between the shock normal 

and V measures the obliqueness of the shock at the point 
sw 

of reflection. The angle V is non zero only if the in¬ 
cident particle has a non-zero thermal velocity component 
parallel to the field. 

Figure A4.1c. Ratio of the particle energy after and be¬ 
fore reflection as a function of shuck obliquity. 

[the above is paraphrased from Sonnerup (1969)] 
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Figure A4.1e shows the energy ratio as a function of shock 
obliquity ( 4 >) for the values - S = V "=0 ^= 0 . 

It should be noted that the singularity in equation A4.8 at 
y>-, 'fl/j. corresponds to the case in which transformation to 
the frame in which the E field is zero is not possible. That 
is, the singularity is artificial, physically, and simply 
means that no finite velocity vector v t (see fig. A4.1) can 
be found that makes the solar wind velocity parallel to *B in 
the moving frame. It can be seen from examining the curve in 
figure A4.I that the Sonnerup model predicts very high values 
of energization on the dusk side of the magnetosphere and rel¬ 
atively small ones on the dawn side. Although electrons of 
energy greater than 40 keV have been observed [ Anderson,(1968)] 
and more recently protons of energy 30 to 100 keV [Lin et al., 
(1974)] in the bow shock regions, no such asymmetry has been 
observed. The dawn side energy increases from the model are 
much too low to account for the observed high energy particles 
in this region. On these bases one must conclude that the 
Sonnerup model is probably not correct. 
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A4.2 Fermi acceleration mechanism-[Jokipii,(1966)] 

The basic idea of the Jokipii model of Fermi accel¬ 
eration at the bow shock is that quasi-thermalized parti¬ 
cles escape from the shock region into a region that is 
bounded on one end by the bow shock and on the other by a 
magnetic irregularity that forms a mirror point upstream 
of the bow shock, (see figure A4.2) Under the assumption 
that the particle velocity parallel to the field is much 
greater than the velocity of the moving irregularity, the 
number of collisions per second that the particle makes with 
the mirror point is 


V / 2 L (t) A4.10 

P ii 


where 

L(t) = L q - Vt A4.ll 

and the gain in momentum on each collision is 

Ap = m V A4.12 

This treatment differs slightly from Jokipii's but the 
differences are primarily notational; for instance in fig- 
gure A4.2, and t 6 =0. 
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Figure A4.2 Jokipii model of the Fermi mechanism for 
accelerating particles at the bow shock. The trapping 
region is defined by the bow shock at one end and a mag¬ 
netic mirror at the other. The magnetic irregularity moves 
toward the shock with velocity c• 
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From equations A4.10 and A4.12, the time 
gain of momentum is 


(t) 


4^ = mV V/2L =MV /2L 

at p „ p 


u = V /V = cos(pitch angle) 
P. P 


p = m V 


but since. 


. 2 4 2 2 . 1/2 
W = (m e + p c ) 


2,2 22 

p=W/c -me 


f -«/' 2 & 


rate of 


A4.13 


A4.14 


A4.15 


A4.16 


A4.17 


A4.18 


A4.19 
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Comparing equations A4.19 and A4.14 we’ see that 

~ = (iUV/L (t) ) w[l- (m^c 4 /W 2 )j A4.20 

which differs from Jokipii's result only in notation. This 

differential equation was derived under the condition that 

V < < V . If the magnetic moment of the particle is con- 
P 

served in collisions with the shock and mirror then as the 
particle gains energy its pitch angle decreases and it will 
eventually be lost to the solar wind. This sets a limit on 
the energy that can be gained in a process like this and that 
limit depends on tlie dimensions of the trapping region and 
the velocity of the irregularity in the field. One gathers 
from this simplistic approach that a particle could gain par¬ 
allel energy only up to about that associated with the bulk 
flow velocity of the solar wind. Energies observed in the up¬ 
stream regions are much higher than this (up to 100 keV protons) 
(Lin et al. (1974)3 which requires some other process to pro- 
due these particles. Jokipii depends on pitch angle scatter¬ 
ing at a fluctuating shock to avoid losing his particles too 
quickly. Lin et al. [1974] have pointed out that this accel¬ 
eration mechanism predicts a hardening of the energy spectrum 
toward the wings of the shock due to the longer time that a 
trapping region ic in contact with the shock which is not ob~ 


;erved. 



169 


Acknowledgements 

I wish to thank my thesis advisor, Dr. John W. Freeman,Jr., 
for giving ire the opportunity to work with a unique project, 
the Suprathermal Ion Detector Experiment. Dr. Freeman's pat¬ 
ience and understanding have been indispensable in the final 
realization of this thesis. 

I would also like to thank Dr. H. Kent Hills for countless 
hours he has spent discussing the subject matter of this thesis 
with me. His immense knowledge of computers and his astonishing 
attention to.detail have been invaluable in the sucessful com¬ 
pletion of this work. 

I have benefitted greatly from my professional association 
and personal friendship with Dr. Richard R. Vondrak who has con¬ 
tributed much of his time to valuable discussions concerning 
this thesis. | 

I would like to thank Tom Rich for his help in computer 
programming and Dave Nystrom and John McGarity for maintaining 
the sanity of the SDS 910 computer during the past three years. 

I have enjoyed knowing and working with the other members 
of the SIDE group, past and present, namely, Dr. Martha Fenner, 
my good friend and ex-office mate, Drs. R.A. Lindeman and 
R.A. Medrano, Mohammed Ibrahim, Dave Hardy, Henry Schneider and 


Dr. Jiirg Me is ter. 





170 


Thanks also go to Steve Bowling and Mike Geis for many 
valuable discussions. 

Finally, I would like to thank my wife Jackie for unself¬ 
ishly spending many nights as a computer widow and for 
typing the final copy of this thesis. 

This research has been supported in part by the 
National Aeronautics and Space Administration under 


contract NAS 9-5911. 




Bibliography 


171 


Anderson, K. A., Energetic Electrons of Terrestrial Origin 
Upstream in the Solar Wind, JGR, 73., 2387, 1968. 


Asbridge, J. R., Bame, S. J., Strong, I. B. 
of Protons from the Earth's Bow Shock, 


, Outward Flow 
JGR, 73 , 5777, 1968. 


Axford, W.I., The Interaction Between the Solar Wind and the 
Earth's Magnetosphere, JGR, 67 , 3791, 1962. 

Bowling, S. B., The Magnetic Field Structure of the Plasma 
Sheet at 30 Re, M.S. Thesis, Rice U., 1973. 


Colburn, D.S., Currie, R.G., Mihalov, J.D., Sonett, C.P., 

Diamagnetic Solar Wind Cavity Discovered Behind the Moon, 
Science , 158 , 1040, 1967. 

Colburn, D.S., Mihalov, J.D., Sonett, C.P., Waid, S.H., The 
Lunar Cavity in the Solar Wind, Trans. Amer. Geophys. 

Union , 49 , 234, 1968. 

Colburn, D.S., Description of Ames Magnetometer Data from 

Explorer 33 and Explorer 35 deposited in the Data Bank, 
NASA/Ames Research Center preprint, 1969. 

Dyal, P., Parkin, C.W., Colburn, D.S., Schubert, G., Lunar 
Surface Magnetometer Experiment, Apollo 16 Preliminary 
Science Report, 11-1, 1972. 

Dyal, P., Parkin, C.W., Sonett, C.P., DuBois, R.L., Simmons,G., 
Lunar Portable Magnetometer Experiment, Apollo 16 Pre : 
liminary Science Report , 12-1, 1972. 

Fairfield, D.H., Bow Shock Associated Waves Observed in the 
Far Upstream Interplanetary Medium, JGR , ,74,_14, 3541, 

1969. 

Fenner, M.A., Magnetosheath Plasma at 60 Re, M.S. Thesis, 

Rice U., 1971. 

Fenner, M.A., Freeman, J.W., Jr., Hills, H.K., The Electric 

Potential of the Lunar Surface, Proceedi ngs of the Fourth 
Lunar Science Conference , 3, 2877, 1973. 

Fenner, M. A., Observations of Magnetosheath Plasma at the 
Lunar Orbit, Ph.D. Thesis, Rice U., 1974. 

Frank, L.A. and Shope, W.L., A Cinematographic Display of 

Observations of Low Energy Proton and Electron Spectra 
in the Terrestrial Magnetosphere and Magnetosheath and 
in the Interplanetary Medium, Trans. Am. G. U. , ji9, _1, 

279, 1968. 

Frank, L.A., On the presence of Low Energy Protons (5^-E^50 kev) 
in the interplanetary medium, JGR , 75./ 707, 1970. 






172 


Freeman, J.W., Jr., The ALSEP Lunar Ionosphere Detector, Rice U. 
preprint, 1969. 


Freeman, J.W., Jr., Balsiger, H., Hills, H.K., Preliminary 
Results from the Lunar Ionosphere Detector, Apollo 12 
Preliminary Science Report , 83, 1970. 


Freeman, 

J.W., 

Jr., Energetic Ion Bursts on the 

Nightside of 

the 

Moon, 

JGR, 77, 

_lr 239, 1972. 


Geiss, J 

., Beuhler, F., 

Cerutti, H., Eberhardt, 

P., Filleux, Ch 


Solar Wind Composition Experiment, Apollo 16 Preliminary 
Science Report , 14-1, 1972. “ 

Goldstein, H., Classical Mechanics , Addison Wesley Pub. Co., 
Reading, Mass., 1950. 

Greenstadt, E.W., Green, I.M., Inouye, G.T., Correlated Magnetic 
Field and Plasma Observations of the Earth's Bow Shock, 

JGR, 73, 1, 51, 1968. 

Heppner, J.P., Sugiura, M., Skillman, T.L., Ledley, B.G., 

Campbell, M., OGO-A Magnetic Field Observations, JGR , 72 , 
21, 5417, 1967. 

Hills, H.K., ALSEP/SIDE Ephemeris Computations, Rice U. pre¬ 
print, 1969. 

Hills, H.K., Meister, J., Vondrak, R.R., Freeman, J.W., Jr., 

The Suprathermal Ion Detector Experiment, Apollo 15 
Preliminary Science Report , 12-1, 1971. 

Hundhausen, A.J., Composition and Dynamics of the Solar Wind 
Plasma, Rev, of Geophys. and Space Physics , J3, _4, 729, 

1970. 

Johnson, F.S., Evans, D.E., Caroll, J.M., Cold Cathode Gauge 
Experiment, Apollo 14 Preliminary Science Report , 185, 

1971. 

Jokipii, J.R., A Model of Fermi Acceleration at Shock Fronts 
with an Application to the Earth's Bow Shcok, Astrophvs . 

J. , 143,_3, 961, 1966. 

Jokipii, J.R., Correlation of ^30 keV Electron Pulses and 

Magnetic Fields in the Magnetosheath and Beyond, JGR , 73 , 

2, 931, 1968. 

Lin, R.P., Meng, C.I., Anderson, K.A., 30 to 100 keV Protons 

Upstream From the Earth's Bow Shock, JGR , 79 , _4./ 4 89 , 1974. 




173 


Lindeman, R.A., Recurring Ion Events at the Lunar Surface, 

M.S. Thesis, Rice U., 1971. 

Lindeman, R.A., Observations of Ions from the Lunar Atmosphere, 
Ph.D. Thesis, Rice U., 1973. 

Lyon, E.P., Bridge, H.S., Binsack, J.H., Explorer 35 Plasma 

Measurements in the Vincinity of the Moon, JGR , 72 1 6113, 
1967. 

Manka, R.H., Lunar Atmosphere and Ionosphere, Ph.D. Thesis, 

Rice U., 1972. 

McKenzie, J.F., Hydromagnetic Wave Interaction with the Mag¬ 
netopause and the Bow Shock, Planet, and Space Scl. , 18 , 

_1, 1-23, 1970. 

Medrano, R.A., Unusual Solar Wind and Solar Proton Events Ob¬ 
served on the Lunar Surface, Ph.D. Thesis, Rice U., 1973. 

Ness, N.F., Scearce, C.S., Seek, J.B., Initial Results of the 
IMP-1 Magnetic Field Experiment, JGR , 69 , 17 , 3531, 1964. 

Ness, N.F., The Earth's Magnetic Tail, JGR , 70 , 2989, 1965. 

Ness, N.F., Behanon, K.W., Cantarano, S.C., Scearce, C.S., 
Observations of the Earth's Magnetic Tail and Neutral 
Sheet at 510,000 Kilometers by Explorer 33, JGR , 72 , 

927, 1967. 

Ness, N.F., The Geomagnetic Tail, Rev, of Geophysics , 7, 1, 2, 

97, 1969. 

O'Brien, B.J., Reasoner, D.L., Charged Particle Lunar Envir¬ 
onment Experiment, Apollo 14 Preliminary Science Report , 

193, 1971. 

Rich, F.J., Plasma in the Geomagnetic Tail Observed by the 

Charged Particle Lunar Environment Experiment, Ph.D. Thesis, 
Rice U., 1973. 

Scarf, F.L. et al.. Direct Correlations of Large-Amplitude 
Waves with Suprathermal Protons in the Upstream Solar 
Wind, JGR , 75 , 34 , 7316, 1970. 

Sonnerup, B.U.O., Acceleration of Particles Reflected at a 
Shock Front, JGR , 74, 5, 1301, 1969. 

Smith, E.J., Holzer, R.E., McLeod, M.G., Russell, C.T., Magnetic 
Noise in the Magnetosheath in the Frequency Range 3-300 Hz, 
JGR, 72, 19, 1967. 



174 


Spreiter, J.R., Summers, A.L., Alksene, A.Y., Hydomagnetic 
Flow around the Magnetosphere, Plan, and Spac. .ScjLj,, 

14, 223, 1966. 

Snyder, C.W., Clay D.R., Neugebauer, M., initial Result from 
AJjSEP Solar Wind Spectrometer (abstract) , Trans. Am. 

Geo. U., EOS , 51 , 407, 1970. 


Snvder C.W., Clay, D.R., Neugebauer, M., Solar Wind Spectrom¬ 
eter Experiment, Apollo 12 Preliminary Sc i ence Reporjt , 

75, 1970. 


Van Allen, J.A., 
of Energetic 
Polar Caps, 


Fennell, J.F., Ness, N.F., Asymmetric Access 
Solar Protons to the Earth's North and South 
JGR, 76, 19, 4262, 1971. 


Wolfe, J.H., McKibbin, D.D., Pioneer 6 Observations of a 

Steady State Magnetosheath, Plan, and Spac. Sci., 1°., 
1968. 


953, 


Wolfe, J.H., Intriligator, D.S., 
with the Geomagnetic Field, 
511-596, 1970. 


The Solar Wind Interaction 
Spac. Sci. Rev. , 10, 4, 


Zhigulev, V.N., Romeshevskii, E.A., Concerning th f J n t®raction 
y of Currents Flowing in a Conducting Medium with the Earth s 
Magnetic Field, Soviet Phv. Dokl. , 4859, 1959. 



